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Abstract 



Positron and Electron Cosmic Rays represent just a fraction of the Cosmic Ray species 
that arrive to the Earth. These particles have the potentiality to only reveal nearby sources 
because they are affected more by energy loss processes than other types of cosmic-rays. 

Galactic positrons are mainly produced by the interaction of nuclei cosmic rays with 
the interstellar gas. Cosmic rays observations in the positron-electron signal reveal the 
possible presence of an unexpected feature for energies above 10 GeV: the positron fraction 
appears larger in the high-energy range with respect to current theoretical predictions. 
During the years, many explanations have been proposed to elucidate this feature. Those 
are based on new physics, such as dark matter particles annihilations, or aditional standard 
astrophysics processes. 

Cosmic-ray experiments as HEAT, AMS, PAMELA among others, have provided 
high-statistic positron-electron data for range of energies from ^ 500 MeV to 100 GeV. That 
gives important restrictions in the theoretical constraints for positron-electron cosmic-ray 
physics. The contribution of dark-matter to the signal, which would appear as deviations 
with respect to the standard prediction, and some of the dark-matter properties would 
be also constrained. A detailed study of the positron background is important to esti- 
mate the capabilities to discriminate a possible new signal present in the experimental data. 

The propagation of cosmic-ray in the galactic environment can be treated in many 
ways. In the case of positrons and electrons, we model the propagation according to the 
Two Zone Propagation Model in which sources related to dark-matter annihilation and 
secondary production have been considered. The positron electron transport equation 
is solved by analytical methods taking into account the uncertainties related to the 
propagation. In addition, the energy spectra of positrons and electrons are described 
by nuclear and particle physics, where we study different production models and their 
uncertainties. 



ABSTRACT 



The feature in the positron fraction, originally seen by the HEAT experiment, can 
be reproduced in the dark-matter annihilation scenario. As well, we showed the effect of 
propagation uncertainties on the dark-matter signal and we studied the potentiality to 
discover a new signal for AMS02 and PAMELA experiment. The secondary production of 
positrons was also studied, taking care on the theoretical uncertainties related to nuclear 
cross section and propagation. The estimations for the positron flux reproduce current 
available experimental data. The positron fraction is calculated on the basis of our results 
on positron flux and fits performed on electron flux data. We obtain that, depending on 
the electron flux used, the fraction may sizeably change in the high energy range stressing 
more or less the necessity of a "positron excess" feature. 

Finally, we give promising results to disentangle a dark-matter component signal from 
the positron background by studying the uncertainties in the positron electron propaga- 
tion. As well, from the study of secondary positrons, we reproduced the observations and 
stressed the importance of the electron signal. Furthermore, PAMELA observations and 
the forthcoming AMS02 mission will soon allow much better constraints on the cosmic-ray 
transport parameters, and are likely to drastically reduce theoretical uncertainties. 
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Sintesi 



I raggi cosmici di elettroni e di positroni costituiscono solo una frazionc dci raggi cosmici 
che raggiungono la Terra. Queste particelle sono potenzialmente in grado di rivelare le 
sorgenti vicine perche sono piu affette da process! di perdita di energia rispetto ad altri 
tipi di raggi cosmici. 

I positroni galattici sono principalmente prodotti da interazioni dei raggi cosmici 
nucleare con il gas interstellare. Le osservazioni di positroni ed elettroni provenienti da 
raggi cosmici evidenziano una possibilc pcculiarita per energie superiori ai 10 GeV: la 
frazione di positroni (positron fraction) e maggiore rispetto all'attuale predizione teorica in 
questo intervallo di energie. Nel corso degli anni, molte spiegazioni sono state proposte per 
chiarire questa particolarita. Questo e alia base di una nuova fisica, come I'annicliilazione 
della materia oscura, o ulteriori processi astrofisici. 

Gli esperimenti sui raggi cosmici, tra cui HEAT, AMS e PAMELA, hanno fornito grandi 
quantita di dati su elettroni e positroni nell'intervallo di energie che va da ~ 500 MeV a 
100 GeV. Questo impone importanti restrizioni sulla teoria dei raggi cosmici. II contributo 
da materia oscura a questo segnale, che appare come una deviazione rispetto alia previsione 
standard, ed alcune proprieta della materia oscura potrebbero anche essere limitati. 
Un dettagliato studio del fondo di positroni e importante per stimare la capacita di 
discriminare un possibile nuovo segnale nei dati sperimentali. 

La propagazione dei raggi cosmici galattici nell'ambiente puo essere trattata in molti 
modi. Ncl caso di elettroni c positroni, abbiamo modcUizzato la propagazione secondo il 
ModcUo di Propagazione a Due -Zone { Two Zone Propagation Model), alPintcrno del quale 
sono state considerate le sorgenti relative all'annichilazione ed alia produzione secondaria. 
L'equazione per il trasporto di elettroni e positroni e stata risolta con metodi analitici, 
tenendo conto delle incertezze relative alia propagazione. Inoltre, gli spettri in energia di 
elettroni e positroni sono stati descritti dalla fisica nucleare e da quella delle particelle, 
dove sono stati studiati diversi modelli di produzione e le incertezze di questi. 
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La peculiarita nella frazione di positroni, originariamente osservata neU'esperimento 
HEAT, puo essere riprodotta nel contesto di annichilazione della materia oscura. Inoltre, 
abbiamo mostrato gli effetti delle incertezze nella propagazione del segnale a causa di 
materia oscura cd abbiamo studiato la potenzialita di rivelare un nuovo segnale ncgli esper- 
imcnti PAMELA cd AMS02. La produzione secondaria di positroni e stata anchc studiata, 
in considcrazionc dcUc incertezze relative alia sezioni d'urto tcorichc c di propagazione. 
Le stime per il flusso di positroni riproduce i dati sperimentali disponibili. La frazione di 
positroni e stata calcolata sulla base sia dei nostri risultati nel calcolo del flusso di queste 
particelle che dei fit dal flusso di dati di elettroni. In conclusione, abbiamo mostrato che, 
a seconda del flusso di elettroni utilizzati, la frazione puo variare in modo signiflcativo 
nell'intervallo di alta energia accentuando una maggiore o minore necessita di un "eccesso 
di positroni". 

Infine, abbiamo ottenuto risultati promettenti per distinguere un componente del 
segnale dovuta a materia oscura, in considerazione delle incertezze del fondo nella 
propagazione. A sua volta, dallo studio di positroni abbiamo riprodotto le osservazioni e 
abbiamo sottolineato I'importanza del segnale degli elettroni. Le osservazioni di PAMELA 
e la futura missione AMS02 permettera delle migliori restrizioni nei parametri di trasporto 
dei raggi cosmici, che permetteranno di ridurre le incertezze teoriche. 
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Resume 



Les rayons cosmiques do positrons ct d'clcctrons nc constituent q'unc fraction dcs rayons 
cosmiques qui arrivent a la Terre. Ces particules ont le potentiel de seulement reveler des 
sources proches parce qu'elles sont affectees par des processus de perte d'energie plus que 
d'autres classes de rayons cosmiques. 

Les positrons galactiques sont principalement produits par I'interaction de rayons 
cosmiques nucleaires avec le gaz interstellaire. Les observations des rayons cosmiques dans 
le signal positron-electron montrent la possible presence d'unc particularite inattendue 
aux energies au dessus de 10 GeV. La fraction dcs positrons {positron fraction) semble 
plus haute a hautes energies que les actuelles predictions theoriques indiquent. Pendant 
des annees, beaucoup d'explications ont ete proposees pour elucider cette particularite. 
Celles-ci sont basees sur de la physique nouvelle, comme I'annihilation de matiere foncee, 
ou sur des processus astrophysiques additionnels. 

Des experiences de rayons cosmiques comme entre autres HEAT, AMS, PAMELA, ont 
offert une grande quantite de donnees sur les positrons-electrons dans I'intervalle entre 
~ 500 MeV et 100 GeV. Ceci impose de grandes restrictions aux limites theoriques pour 
la physique de rayons cosmiques. La contribution de la matiere foncee a ce signal-ci, qui 
apparaitrait comme une deviation de la prevision normale, et quelques proprietes de la 
matiere foncee, seraient aussi restreintes. Une etude detaillee du fond de positrons est 
importante pour estimer les capacites pour discriminer un possible nouveau signal existant 
dans les donnees experiment ales. 

La propagation do rayons cosmiques dans Tatmosphcrc galactiquc pout otro traito do 
beaucoup de manieres. Dans le cas des positrons et des electrons, nous avons modele 
la propagation en accord avec le Modele de Propagation a Deux Zones {Two-Zone 
Propagation Model) dans lequel les sources concernant I'annihilation et la production 
secondaire ont ete considerees. L'equation de transport des positrons-electrons a ete 
resolue par des methodes analytiques en prenant en consideration les incertitudes del la 
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propagation. Aditionnellement, les spectres d'energie des positrons et des electrons ont ete 
decrits par la physique nucleaire et la physique des particules, dans les quelles on a etudie 
de differents modeles de production et des incertitudes de ceux-ci. 

La particularite dans la fraction de positrons, initiellement vue dans I'cxpcrience 
HEAT, pcut etre reproduite dans le scenario dc I'annihilation dc maticrc fonccc. En outre, 
nous avons montrc les cffcts des incertitudes dc propagation du signal, ctant donnc la 
matiere foncee et nous avons etudie la possibilite de decouvrir un nouveau signal dans les 
experiences AMS02 et PAMELA. Le potentiel secondaire de positrons a aussi ete etudiee, 
considerant des incertitudes theoriques relatees aux sections efficaces et a la propaga- 
tion. Les estimations pour le flux de positrons reproduisent les donnees experimentales 
disponibles. La fraction dc positrons est calculcc sur base dc nos rcsultats dans le calcul 
des flux dc cc\ix-ci ct des ajustcmcnts faits sur les donnccs dc flux d'clcctrons. Nous avons 
obtenu que, suivant le flux d'clcctrons utilise, la fraction pcut considerablcnicnt changer 
dans I'intervalle des hautes energies accentuant plus ou moins la particularite d'un "exces 
de positrons". 

Finalement, nous donnons des resultats prometteurs pour decouvrir un composant du 
signal, etant donne matiere foncee, du fond en etudiant des incertitudes de la propagation. A 
son tour, a partir de I'etude de positrons secondaires, nous avons reproduit les observations 
et nous avons insistc sur I'importancc du signal d'clcctrons. En outre, les observations de 
PAMELA et de la future mission AMS02 permettront d'obtenir de meilleures restrictions 
aux parametres de transport des rayons cosmiques, et seront prometteuses pour reduire 
des incertitudes theoriques. 
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Resumen 



Los rayos cosmicos do positroncs y clcctroncs constituycn solo una fraccion dc los rayos 
cosmicos que Uegan a la Tierra. Estas particulas tienen la potencialidad de revelar solo 
fuentes cercanas porque son afectadas en mayor medida por procesos de perdida de energia 
que otras clases de rayos cosmicos . 

Positrones galacticos son principalmente producidos por la interaccion de rayos 
cosmicos nucleares con el gas interestelar. Las observaciones de rayos cosmicos en la senal 
positron electron mucstran la posiblc presencia de una inesperada peculiaridad en cl rango 
por sobre de los 10 GeV: la fraccion de positrones {positron fraction) se muestra mas grande 
en el rango de altas energias que las actuales prediciones teoricas. Durante anos, muchas 
explicaciones ban sido propuestas para dilucidar esta peculiaridad. Estos estan basados en 
nueva fisica, como la aniquilacion de materia oscura, o en procesos astrofisicos adicionales. 

Experimentos de rayos cosmicos como HEAT, AMS, PAMELA entre otros, ban 
brindado gran volumen de datos de positrones y electrones en el rango que va desde 
^ 500 MeV hasta los 100 GeV. Esto impone grandes restricciones a los li'mites teoricos 
para la fisica de rayos cosmicos. La contribucion de la materia oscura a esta serial, que 
apareceria como desviaciones respecto a la prediccion estandar, y algunas propiedades de 
la materia oscura seri'an tambien restringidas. Un estudio detallado del fondo de positrones 
resulta importante para estimar las capacidades de discriminar una posible nueva senal 
existente en los datos experimentales. 

La propagacion de rayos cosmicos en el ambiente galactico puedc scr tratado de 
muchas formas. En el caso dc positrones y electrones, hemos modelado la propagacion de 
acuerdo al Modelo de Propagacion a Dos Zonas {Two-Zone Propagation Model) dentro 
del cual las fuentes relacionadas a la aniquilacion y a la produccion secundaria ban sido 
consideradas. La ecuacion de transporte de positrones y electrones ha sido resuelta por 
metodos analiticos tomando en cuenta las incertezas relacionadas a la propagacion. En 
anadidura, los espectros de energia de positrones y electrones ban sido descritos por fisica 



vii 



RESUMEN 



nuclear y de parti'culas, donde han sido estudiados distintos modelos de produccion y las 
incertezas de estos. 

La peculiaridad en la fraccion de positrones, originalmente observados en el experi- 
mento HEAT, puede ser reproducida en el escenario de la aniquilacion de materia oscura. 
Ademas, hemos mostrado los efectos de las incertezas en la propagacion dc la scfial debida 
a materia oscura y hcmos cstudiado la potcncialidad dc dcscubrir ima nucva serial cn los 
experimentos AMS02 y PAMELA. La produccion secundaria de positrones tambien ha 
sido estudiada, considerando las incertezas teoricas relacionadas a las secciones eficaces 
y a la propagacion. Las estimaciones para el flujo de positrones reproducen los datos 
experimentales disponibles. La fracion de positrones es calculada en base a nuestros 
resultados en el calculo de los flujos dc cstos y dc ajustcs hechos sobrc datos de flujos de 
electrones. Obtuvimos que, dcpcndicndo del flujo dc clcctrones usado, la fraccion pucde 
cambiar considcrablcmcntc cn cl rango dc altas cncrgias acentuando en mayor o menor 
medida la peculiaridad de un "exceso de positrones". 

Finalmente, damos prometedores resultados para desentranar una componente en la 
senal debida a materia oscura del fondo, estudiando las incertezas en la propagacion. 
A su vez, del estudio de positrones secundarios, hemos reproducido las observaciones y 
hecho hincapie en la importancia de la serial de electrones. Ademas, las observaciones de 
PAMELA y dc la futura mision AMS02 pcrmitiran mejores restriciones a los parametros 
de transporte de rayos cosmicos, y ser an prometedoras para reducir las incertezas teoricas. 
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Chapter 1 



Introduction 

In the last decades, the astrophysical and cosmological evidence of Dark 
Matter and Dark Energy have created a revolution in the field of fundamental 
physics. Some models, which are candidates to replace the current Standard 
Model of particles physics, predict dark-matter particle candidates. In similar 
way, the cosmic-ray physics have been stimulated with observations that 
show a much active universe that was thought. As well, cosmic-ray are 
promising proves in the understanding of the local environment. Those are 
genuine samples of the matter composition of the galaxy. In addition, the 
antimatter cosmic-rays component, which are less abundant than matter 
cosmic-rays component, gives crucial clues regarding to the non-standard 
contribution to the cosmic-ray signal. 

1.1 Dark Matter 

The definition of Dark Matter (DM) comes out from the fact that this kind of matter 
does not emit or absorb electromagnetic radiation at any wavelength, its gravitational 
interactions dominate on scales from tiny dwarf galaxies, to large spirals such as the Milky 
Way, to clusters of galaxies, to the largest scales until now observed. 

Spiral galaxies support the hypothesis of DM in which star dynamics suggest the 
presence of additional mass which is not detectable using electromagnetic radiation 
detection. Moving to larger scales, such as galaxy clusters, the evidence of DM comes from 
different experimental methods, such as gravitational lensing, X-ray gas temperatures and 
the motion of cluster member galaxies. In general, depending on the scale to which we look 
at, different methods of measuring directly or indirectly the presence of DM are employed 
(Figure OJ. 
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Figure 1.1: Left panel: Indirect evidences of Dark Matter present at all scale. Right panel: Pie 
chart of composition of the observable Universe, Dark Matter is the 23% at current 
time [3]. 



As well, the Dark Matter component in the Standard Cosmological Model is fundamen- 
tal to explain many of the current observations. The combined WMAP, SNIa and galaxy 
clusters measurements suggest the presence of DM on cosmological scales. Its contribution 
to the Universe energy content is [U [2] 

r!DM^0.23, (1.1) 

which corresponds to the 23% of the total energy density. A satisfactory description of 
most cosmological observations is obtained by the so called "ACDM model", which comes 
out from the best fit of the combined data analysis. 



1.1.1 Dark Matter particle candidates 

Potentially the only indication compatible with cosmological measurements is that dark 
matter is composed of non-baryonic, neutral and weakly interacting particles. In the liter- 
ature several candidates were proposed, the most relevant are: 

Standard Model Neutrinos: The existence of a relic sea of neutrinos in number only 
slightly below that of relic photons that constitute the CMB, is a generic prediction of the 
standard hot Big Bang model. Their contribution to the matter density of the Universe is: 

= (1.2) 

The requirement that Q^, < 57^ ~ 0.3 imposes stringent limits on their masses. Indeed 
dark matter particles with a large velocity dispersion such as that of the neutrinos affect 
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the evolution of the cosmological perturbations. This leads to a top-down scenarios which 
is not supported by the present observations since the galaxies seems older than cluster. 

Heavy Neutrinos: The allowed mass range is bounded from below from the Lee- 
Weinberg limit |4j, which was > 2 GeV at the time of their work. The current limits 
have been updated to > 1.3 — 4.2 GeV for Dirac neutrinos and > 4.9 — 13 GeV for 
Majorana ones [5]. A more stringent bound comes from colliders: neutrinos lighter than 
45 GeV are excluded by the total decay width of the Z boson. For very heavy neutrinos 
the Yukawa coupling to the Higgs boson would be so strong that perturbative calculation 
become non- reliable [6,^ and partial waves unitarity has to be imposed [J, leading to a 
stringent mass upper bound, < 3 TeV. In the allowed mass range the cosmological 
properties can be interesting, however their interactions are quite strong, therefore they 
are mainly excluded by direct detection bounds and lead to a low relic abundance. 

Sterile Neutrinos: Those are similar to Standard Model neutrinos, but without Stan- 
dard Model interaction apart from mixing [8J. Stringent cosmological and astrophysical 
constraints on the sterile neutrinos come from the analysis of their cosmological abundance 
and the study of their decay products |9j . 

Axions: Originally those have been introduced to explain the so called strong CP 
violation problem |101 lll|. From different searches, it is expected that axions are very light 
and extremely weakly interacting with ordinary particles, which implies that they were 
not in thermal equilibrium in the early Universe. The calculation of the relic density is 
uncertain, nevertheless it is possible to find some range where axions satisfy the present 
constraints and represent a possible DM candidate |l2j . 

Supersymmetric Particles: In models conserving i?-parity the lightest supersymmetric 
particle (LSP), such as the neutralino, sneutrino, gravitino or axino could provide the right 
amount of dark matter density in the Universe. Which particle is the LSP, it depends on 
the supersymmetric models and on how supersymmetry is broken. 

Kaluza-Klein states: Those are excitations of the Standard Model fields which appears 
in models with extra dimensions [T31 HH HH [Hj . 
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Super— heavy dark matter: It is composed by heavy stable particles with a mass in the 
range of 10^^ to 10^^ GeV. These particles leads to scenarios for production of nonthermal 
dark matter [T7] . 

Light scalar dark matter: A class of fermionic dark matter candidate, Lee and 
Weinberg concludes that relic density arguments preclude such a WIMP with a masses 
less than a few GeV [T5]. The dark matter is form by other stable dark matter species, as 
in the case of SUSY N = 2 theories, where both types coexist at the same time. 

Dark matter from Little Higgs models: These models have been proposed in order to 
stabilize the weak scale and to solve the hierarchy problem as an extension of the Standard 
Model. This class of models possess discrete symmetries which results in the existence of 
stable weakly interacting particles [IHl [20] . 

1.2 Cosmic Rays 

The cosmic rays (CR) are generically described as charged particles that travel across 
the Universe. Their origin is usually associated to the nuclear activity present in stars, 
galaxies, etc. Depending on the processes in which they are involved, they are present at 
different energies scales usually going from the eV- up to PeV-scale. 

A first classification of CR regards their production location: 

• Solar CR. 

Also known as solar energetic particles (SEP), these are cosmic rays that originate 
from the Sun. The average composition is similar to that of the Sun itself [H]. 

• Galactic CR. 

This type consists of those cosmic rays that enter the solar system from the outside. 
They are high-energy charged particles composed of protons, electrons, and fully 
ionized nuclei of light elements. 

• Extragalactic CR. 

Unlike solar or galactic cosmic rays, little is known about the origins of extragalactic 
cosmic rays. This is largely due to a lack of statistics: only about 1 extragalactic 
cosmic ray particle per square meter per year reaches the Earth's surface. 
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Many other ways to classify them are possible. The nature of the cosmic-rays also 
establish a good rule to classify them. Generally, those are separated into: Nuclei CR and 
Electron CR. Each category involves different features. Nuclei CR are less affected by energy 
losses, unlike electron CR that can cover shorter distances. This characteristic allows to 
nuclei CR to travel longer distances and increase their chances to interact with the medium. 

On the other hand, we have the antimatter CR, which are less abundant than matter 
CR. The origin of this type of particles are not well understood, even though a fraction of 
those are related to spallation process between CR and the interstellar gas. 

1.2.1 Electron and positron cosmic— rays. 

Let us focus on the electron and positron CR species. As we said before, those particles may 
travel shorter distance in opposition to nuclei ones. This would give detailed information 
about the Earth's local environment and/or the presence of exotic component, because a 
shorter distance reduces the chances to interact more with the medium. 

The production of this species is generally related to: 

• Supernovae and gas: 

Cosmic-rays are injected into the medium when supernova explosions occur. The 
effect of Shockwaves allows the particles in the medium to gain energy and travel 
across space. 

• Secondary production: 

Primary cosmic-rays interact with the gas present in the interstellar and intergalactic 
space. The spallation processes allow the production of new CR, which start to 
propagate and contribute to the cosmic-ray signal [22] . 

• Pulsars: 

A number of studies have been done involving pulsars as sources of electron 
CR and recently of positrons CR as well ^5] . 

• Exotic sources: 

Dark-matter particle annihilation [T31 and evaporation of Primordial Black 

Holes make possible contributions to the CR signal that come from the 
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Figure 1.2: Electron (left) and positron (right) fluxes versus energy. Some of the lastest exper- 
imental results and the Strong and Moskalenko prediction [391 [40] are shown. The 
data correspond to HEAT [3D], AMS [SUET], CAPRICE [3T], MASS [35], PPB- 
BETS [33 [34] and ATIC-2 [35] experiments. 

galactic halo. 

1.2.2 Positron and electron observation. 

Most of the observations of galactic electron and positron cosmic-rays are performed by 
balloon and space-borne experiments. Some of the balloon experiments are HEAT |30j . 
CAPRICE [3T|, MASS [32], PPB-BETS [31 [M] and ATIC-2 [35 . The most well known 
space-borne experiments are AMS [SS] [37] , which was attached to the space shuttle in one 
of its missions, and PAMELA (35] that is currently on flight since 2006. 

The current status of electron and positron flux observations is shown in Figure 11.21 
We observe that electrons arrive more abundant than positrons. As well, the effect of 
Solar activity on the low-energy range affect the fluxes reducing and deforming them 
with respect to the shape at the Solar System boundary. The theoretical predictions 
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done by Strong and Moskalenko [39 show agreement with respect to the observation 
in both signal. A different situation occurs when the positron fraction, the ratio be- 
tween the positrons and the total amount of positron and electrons, is plotted. In 
Figure II .31 the positron fraction for different experiment is reported. We can observe 
that the Strong and Moskalenko prediction deviates in the high-energy tail. This has 
produced a big revolution in the field, because there are not certain answers to this problem. 

In this thesis, we study the positron and electron cosmic-ray signal. The work is 
oriented to solved and understand the "positron excess" feature. In Chapter [2l we study the 
different production mechanisms, that occur in the galactic environment, related to nuclear 
and particle physics. We study the production of positron and electron in proton-proton 
interactions and from dark-matter annihilation-like processes. In Chapter [sj we review 
the propagation model for ultrarelativistic cosmic-ray and improve the current solutions 
regarding the positron and electron case. 

Chapter |4] is devoted to study the hypothesis that DM annihilation is the responsable 
of the "positron excess" feature. We analyze the DM annihilation signal from the context 
of generic-DM candidates and propagation uncertainties. As generic-DM candidate, 
we consider some DM annihilation channels which correspond to typical signatures of 
candidates proposed by Beyond the Standard Model theories, as SUSY or Kaluza-Klein 
DM particles. The positron background is studied in Chapter [S] We study the secondary 
positrons produced by the interaction of proton and alpha particle CR with the interstellar 
gas present in the galactic plane. As well, we reanalyze the "positron excess" problem using 
the available experimental data. In this case, we found that the positron fraction is sensible 
to small variation in the positron flux. 
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Figure 1.3: Positron fraction versus energy. The current observational status of the positron 
fraction and the prediction done by Strong and Moskalenko [391 I4U] are shown. 
The prediction shows the presence of a "positron excess" in the high-energy tail. 
The data correspond to HEAT [30], AMS [MIEZ], CAPRICE [311, MASS [32] and 
PAMELA [H] experiments. 
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Chapter 2 

Positron and Electron 
production 

The production of electrons and positrons in galactic environment takes 
place in different ways. These processes are mainly described by standard 
particle physics. In this chapter, production processes related to proton- 
proton interactions and Dark Matter annihilation are reviewed and described. 

2.1 Overview 

The production of positrons and electrons is possible through different processes and in 
a wide range of energies. In the observable Universe, electrons are present in every place 
where atoms are. Electrons can be removed from atoms and be expelled into the interstellar 
space as result of supernova explosions, gas ionization processes and pulsars interactions 
with the medium. As well, electrons are produced by CR interaction with the Interstellar 
Medium (ISM) 

The case of positrons is quite different. The principal contribution comes basically from 
CR interaction with ISM |39j . The matter /antimatter abundance asymmetry makes that 
mechanisms similar to the electron's ones become highly suppressed. 

In the range of energies below 1 GeV, CR measurements, specially electron and positron 
CR fluxes, can roughly be explained with our actual knowledge in nuclear and particle 
physics and nearby sources. Nevertheless, measurements above this limit may suggest 
the presence of undiscovered sources or new physical processes related to production 
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An undiscovered category of sources may lay on Beyond the Standard Model (BSM) 
theories context [12] ■ These theories suggest that exotic CR signals may be related to CDM 
particles, where positrons and electrons are a consequence of their annihilations. In this 
category, sources are described and modeled in terms of initial and final states, formed by 
BSM particles and SM particles respectively. 

In this chapter, electron and positron production related to proton-proton and DM-like 
interaction are reviewed and explained. Let's warn that topics described in this chapter 
are crutial to understand the following ones. 

2.2 Production in Proton— Proton collision 

Physics behind cosmic rays sources are closely related to standard particle physics. In 
fact, most of processes where CR production takes place are explained by the last one. 
Collisions between nuclei CR and ISM are sources of secondary CR. Secondary positrons 
and electrons are also produced in that way. 

The ISM is composed mainly by Hydrogen and Helium with densities of 
nn = 0.9 [part/cm^] and riHc = 0.1 [part/cm^] respectively [i^ liH B51 liS] . 

Nuclei CR are mainly constituted by protons and alpha particles. Many experiments are 
devoted to measure CR fluxes. Generically, interstellar (IS) CR fluxes are parameterized 
following a power-law form [47 : 



which depends on the CR rigidity (TZ = — — ) and Lorentz boost factor (3. 

Ze 

The values of parameters depend on the CR type. In the case of protons, those are: 




(2.1) 



A= (1.94 ±0.13) X 10-*, pi 



0.70 ± 0.52 and pa = 2.76 ± 0.03. 



(2.2) 



The values for alpha particles are: 



A = (7.10 ±0.56) X 10^ pi 



0.50 ± 0.31 and pa = 2.78 ± 0.03. 



(2.3) 
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Positrons (electrons) produced in CR and ISM interactions will depend directly on 
the energy dependence of proton and alpha fluxes. As well, there is an intrinsic positron 
(electron) energy distribution related to the collision itself. The intrinsic distribution is 
out of the astrophysics context, and it can be explained by using just particle and nuclear 
physics. 



2.2.1 Production through mesons decay 

The number of particles produced in hadronic collision, like in the proton-proton case, 
grows faster than in leptonic ones. In general, lighter mesons are produced in big quantities, 
specially pions and kaons. Those have decay modes which produce many types of particles, 
for example, neutral pions are well known as an efficient source of gamma-rays. As well, 
charged pions and kaons are efficient in production of positrons and electrons. 

Positive and negative pions - or kaons - are produced in equal number due to con- 
servation laws. However, positron and electron energy distributions are not equal. The 
difference is related to polarization effects in the muon production (Appendix |b|), the 
left-right asymmetry present in weak interactions is responsable of this. The result is that 
electrons are more energetically produced than positrons. 

In general terms, positron and electron inclusive productions are not so different in 
form; the main difference lays on the muon's polarization effect. Apart of that, those are 
calculated in same way. The first ingredient is the inclusive cross sections (CS) of charged 
mesons, which are used to guide the positrons (electrons) production (through the decay 
into positrons and electrons). 

In the literature, there are many models to explain mesons production. However, we 
decided to use parameterizations of inclusive CS, in order to obtain more realistic results. 
In the following, we explain the most used ones. 



Badhwar Stephens Golden parameterization 

This was proposed by Badhwar, Stephens and Golden in the late 70's [IHIHS]- Its aim is 
to explain and to get good agreement with accelerators and CR data. Basically, it is a 
parameterization of the invariant CS in p-p collisions. 
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Particle 


7r+ 


TT 


K+ 


K- 


A 


153 


127 


8.85 


9.3 


B 


5.55 


5.3 


4.05 


3.8 


r 


1 


3 






C 






2.5 


8.3 


Ci 


5.3667 


7.0334 






C2 


-3.5 


-4.5 






C3 


0.8334 


1.667 







Table 2.1: Parameters for charged Pion and Kaons inclusive CS according to Badhwar-Stephens- 
Golden parameterization [481149) . The parameter units are such that the inclusive CS 
are in units of mb GeV~^ f Equations 12.41 and 12 . lOp . 



Particle 


Exclusive reaction 


mx (GeV/c^) 


7r+ 




1.878 






2.016 


K+ 


pA°K+ 


2.054 


K- 




2.370 



Table 2.2: Exclusive reactions in p-p collisions and the minimum configuration mass mx 
f Equation 12 . 9[1 for charged pious and kaons [50j . 

The parameterization for production of charged pions is: 

e^-^ = AK^'-il ~ exp pt k^) (mb GeV-^ c^) , (2.4) 

dp'^ 

where pt is the transverse momentum, 

q{pt) ^ {Ci + C2 Pt + p't) ^ and j^y^2^lM±Z!^'. (2.5) 

Values for the parameters are given in Table [2~T1 However, it is necessary to remark that 
this is a parameterization of differential inclusive CS in the laboratory frame, that means: 

7Lab = ^and n-' = ^l+j-^J . (2.6) 



The variable which appears in Equation 12. 51 is formally defined as the ratio among 
momenta in the center of mass system (CMS): 



P, 

12 



^ (2 7) 

CMS ' 
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where 

„CMS _ f^Lah a ^Lab 



pY"''' = TLabbfl"" - pLabS^n , (2-8) 



„CMS , ^CMS _ s + — rrix 



We need to specify that m is the produced meson mass and mx is the mass of minimal 
production configuration, which also gives infomation about thresholds for the energy of 
center of mass (ECM) and meson energy flable [2~2|1 . 

A second parameterization is proposed, this time is for charged kaons and it is slightly 
different respect to the pions case: 

e^-^^A{l-S:fexp{-Bpt) (mb GeV-^ c^) , (2.10) 

where A, B and C are constants with values given in Table 12. 11 Let's emphasize that the 
parameterization fits good observational data. 



Tan Ng parameterization 

This is a well known and used parameterization for the study of inclusive CS's in p-p col- 
lisions. It reproduces and predicts with good precision the low and high energy data [51[l50j . 



The invariant CS is generically written as: 

-"^'^"^"^ -f{xn)exp{~A{xn)pt-B{xM (mb GeV^ c^) , (2.11) 



dp 
where 



^CMS 



is the radial scaling variable, which is expressed in terms of meson energy and maximum 
possible energy for the minimal configuration in the CMS (Table 12. 2p . Another scaling 
variable is 

XT = — — , (2.13) 

1 Xni 



where 

^™ = pCMS ' (2-14) 

max 

13 



CHAPTER 2. POSITRON AND ELECTRON PRODUCTION 



Particle 


7r+ 


TT 


K+ 


K- 




163 


163 


7.33 


7.33 




143 


150 


2.61 


6.97 


dn 


9.54 


13.8 


5.36 


11.8 


O.-i 


0.7 


0.5 


0.7 


0.5 




2.72 


3.38 


2.49 


3.27 




3.53 


4.20 


1.96 


2.27 




3.48 


4.82 


1.03 xlO~^ 


4.32 


O7 


1.88 


1.07 


1.67 


4.89 xl0"2 


as 


0.212 


0.984 


0.263 


6.03 


ag 


1.64 xlO-"* 


0.873 


8.24 xlO^ 


22.6 


aio 


0.779 


0.455 


6.41 


0.894 


an 


4.70 


1.71 


4.72 


2.81 


ai2 


7.82 


2.01 


6.09 


0.375 



Table 2.3: Parameter for charged Pion and Kaons inclusive CS. These sets are for Tan-Ng 
parameterization |51l . The parameter units are such that the inclusive CS are in 
units of mb GeV"^ (Equation l2 . lip . 

it helps to extend the parameterization to very low energy scales, almost arriving to 
threshold energy, with pretty good precision. 

The resting terms are: 

/(xr) = ai exp ( - a2a;R)6(a3 - xr) + (ctoo - ai)(l - xr)"'' (2.15) 
A{x'r) = 05 exp ( - aexn) + aj cxp (agXR) 
B{xB.) = a9exp(-aio(xR + aii))(xR + aii)°i^ 

where Q is just the Heaviside step function. 

For our purposes, we are interested in the production of charged pions and kaons, for 
which, the parameter values are given in Table 12.31 

Stecker model 

The Stecker model was created to explain the pions production at very low energies in 
p-p collisions. The principal mechanism is through the production and decay of A. Also, a 
secondary pion production mechanism is related to fireballs production, i.e. thermal pion 
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production from residual collision energy [551 [53] . 

In the p-p collision, production of A baryons produce charged and neutral pions in 
almost the same proportion. Originally, Stecker's model was used to estimate gamma-rays 
fluxes. However, his ideas are used also for the production of charged pions, and then 
positron and electrons. 

The discussed process is: 

p + p — >p+ (A+ ^ n + 7r+) — >p + n + 7r+ , (2-16) 

where A baryons are considered as resonances, i.e. its effective mass varies following a 
Breit-Wigner distribution, 

where — 1.232 GeV and T = 0.118 GeV are the most recent determinations [2]. Some 
physical considerations reduce the niA range to: 



'A 



TOp . (2.18) 



The first bound corresponds to the lowest value of mass able to produce a physical 
pion and the second bound is related to the minimal configuration, i.e. the production of 
a proton with a A and nothing else. 

The A decay into charged pions may happen just in one way, and so, the pion energy 
distribution will be: 

S{ejr — when A is at rest. 

/^,A(e^,eA) = <( 1 (GeV-i) , (2.19) 



27a/3aP* 



moving A and e™'" < < e™' 



where 



is the pion energy in the A rest frame. In the case of a moving A, pion kinematic limits 
are: 

ef^ = ^Aiel - PapD and er"" ^ lAiel + PapD , (2.21) 
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Figure 2.1: Schematic represetantion of pion production resulting from p-p collision according 
to Stecker's model. 



where 



/3a = 



PA 



(2.22) 



the description of moving A is based on the boosted spectra method (see Appendix |a|) . 

A third point to discuss about is the A production in p-p collisions. In the CMS, the 
energy distribution of A baryons is: 



where 



s — nip + 



2V^ 



(2.23) 



(2.24) 



is solution in the CMSfor 2-body production and an isotropic A production is assumed. 



However, we are interested in the solution in the laboratory frame. Using the boosted 
spectra formalism (Appendix |a|), we found the equivalent energy distribution in the this 
frame: 

1 



/a^Pp(£A, Ep) 



VeAe[£S'",£n (GeV-i), (2.25) 



2 7Lab/3Lab P*A 

also the kinematic limits change giving new bounds for the A energy: 



= 7Lab(el - /3LabPA) and e^'^'^ = 7Lab(eA + Ph^hPl) ■ 



(2.26) 



For a fixed A mass, the pion energy distribution is obtained as a composition of pro- 
cesses: 

pLab ^ ™ \ / rhah 



f^^p';{e^,ep,mA)^ l_ d^A /^"pVeA, £p) A,A(e-, ^a) , (2.27) 
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where after integration, we obtain: 
where 



7A(e.) - max (^7^^ '"'"^2^"^ J ' (2-29) 
which is a typical behavior related with the change of reference frame (Appendix |a|) . 

The next step is to include the A-resonance. For that, an average respect to mA is 
needed: 



f^fpie^, ep)^N dmATimA) fi^fpie., Sp, mA) (2.30) 

J„jmi„ 

where A'^ is a normalization constant for J^(mA), in which, 

= - arctan ( ^ ^ j - arctan ( ^ ^ j j . (2.31) 



Unfortunately, this integral cannot be solved analytically and numerical integration 
should be performed [52l |53] . 



Calculation of inclusive cross sections 

Inclusive OS's for positron and electron are calculated as convolutions between mesons 
OS's and their energy decay distributions into positrons and electrons. 

The final positron (electron) inclusive CS is obtained from contributions of pions and 
kaons: 

^ djjja)^ ^ djjja)^ ^^^^^^ 

dSe dSe dSe 

where each contribution is generically calculated as: 

^^^^^^(£e,ep) = J dsx^^^^^^^{sx,Sp) X fe,x{se,ex) (mb GeV-1) , (2.33) 

where X denotes pions or kaons as intermediate particles. Note that meson CS are re- 
lated to the parameterizations previously seen, through an integration over the transversal 
momentum: 



djycr) 
de 



{e,ep) = nj^ (^^) (Pt^Pf^p) ^ (-^ GeV-^) . (2.34) 



P\\ 
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The decay distributions fe^x have been analyticaUy calculated including muon polar- 
ization effects. Pion decays are the simplest ones because they have just one dominant 
decay mode (details in Appendix [c|) . On the contrary, Kaon decays are more complex to 
treat (details in Appendix |d|) . 



2.2.2 Kamae et al. parameterization 

A parameterization for production of positrons, electrons and other particles are proposed 
by Kamae et al. [Ml [55l [56]. The objective of this parameterization is to give an easy way 
to compute and estimate CR fluxes, that comes from ISM interactions with nuclei CR. 

New processes are included, contributions from A(1238) and several hadron resonances 
around 1600 [MeV/c^] make it more accurate in the very low proton energy range. Another 
included process comes from diffractive dissociation which contributes in an intermediate 
energy range. 

For the positron and electron cases, the inclusive differential CS is described by: 

rf(w)ND(£e) ^ ^^y^ ^ p^^^^^ ^ F^u.M (mb) , (2.35) 

alog(£e) 

where y — logio(Tp/lTeV) and x — logiQ(ee/lGeV). 

The function r{y) is a rescaling of the Non-Diffractive (ND) contribution in order to 
reproduce experimental data. -Fnd is the parameterization of ND CS: 

^ND(a;) = flo cxp |-ai(a; - as + 02(0; - 03)^)^1 + (2.36) 
04 exp |— 05(0; — as + ae{x — ag)"^ + 07(2: — ag)^)^! , 

The parameter values and r{y) functions are given in Table 12.61 for positrons and in 
Table [2~5] for electrons. -Fnd.m represents the conservation of energy- momentum: 

-Fnd ki(a;) = p n ^ r Ti • (2-37) 

exp jWND.U-^min - a;)) + 1 exp |WND,h(a; - imaxj) + 1 

where values of parameters are give in Table 12.41 
Diffractive processes are parameterized as follows: 

Fo„(.) . 6. exp + »p { -H ( ij;,"!,^) ) (2.38) 

18 



2.2. PRODUCTION IN PROTON PROTON COLLISION 



Particle 






e~ 
e+ 


-2.6 0.96 X logio {Tp/lGeV} 
-2.6 0.94 X logio {Tp/lGeV} 


20 45 
15 47 



Table 2.4: Kinematics limits for Kamae et al. parameterization for electron and positron pro- 
duction. 



and the resonance contributions are 



X - C2 



(2.39) 



1 + C3{X - C2) + Ci{x - C2Y , 

As before, the energy-momentum conservation is a separate function and is described 



by: 



Fvxix) 



(2.40) 



exp { (Wdiff (X - Lmax)} + 1 ' 

where Wdiff = 75 and L^ax — log^o {Tp/lGeV}. The inclusive CS is composed as before 
rf(w>ND(e( 



. ^ ^ ^Diff/rcs(a;) X Fki(a;) (mb), 
dlog(£e) 

although the rescaling r(y) function is not included here. 



(2.41) 



Kamae et al. parameterizations work well for positron and electrons. However, there 
were changes in the parameter values respect to published ones [ST]. Let's clarify that the 
new values are given in all the tables. 



2.2.3 Production uncertainties 

Since each inclusive CS parameterization is based on physical assumptions and it repro- 
duces experimental data, there are uncertainties related to parameter determination, which 
modify the asymptotical behavior at low and high energy ranges. 



As it can be observed in Figure [2~2l at different proton energies, positron and electron 
CS for the three parameterizations are closely similar in behavior. However, there are 
variations up to 80% at proton energies of 20 GeV, as in the case of Kamae et al. versus 
Tan and Ng parameterizations. Another feature is that Kamae's parameterization estimate 
a smaller electron CS respect to the other two parameterizations. 
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Parameterization for Electrons 




Parameters 


Formulae as fmictions of the proton kinetic energy {y = logjo(rp)) in TeV. 




0,0 


-0.018639(2/ + 3.3) + 2.4315(2/ + 3.3)^ - 0.57719(y + 3.3)^ + 0.063435(y + 3.3)* 




Ol 


7.1827 X 10"*^ - 3.5067 x IQ-'^y + 1.3264 x lO^^^y^ - 3.3481 x lO'V + ■ ■ • 




. . . + 2.3551 X 10-^2/* + 3.4297 x W'^y^ 




02 


r , 2.8924 X iC 
563.91 -362.181ogio {2.7187(2/ + 3.4)} - ^^^^ ^^3^^, 


o 


03 


0.52684 + 0.577172/ + 0.00453367/2 - 0.00890662/=* 


M 


04 


0.36108(?/ + 3.32) + 1.6963(?/ + 3.32)^ - 0.074456(2/ + 3.32)^ - . . . 


3 


... - 0.071455(2/ + 3.32)'' + 0.010473(y + 3.32)^ 






9.7387 X 10-5 + 7.8573 X lO-^logjo {0.0036055(2/ + 4.3)} + .. . 




ae 


0.00024660 , 
y + 4.9390 " 
-273.00 - 106.22 logn, {0.34100(2/ + 3.4)} + 89.0372/ - 12.5462/^ 




07 

ttS 


432.53 -883.991ogio{0.19737(y + 3.9)} - ^y^g^^^g", 
-0.12756 + 0.434782/ - 0.00277972/^ - 0.00830742/^ 






Q (^Q J ^r\a f y ~^ 3.26 ^1 '7c2 rp 1 c rf^ir 
o.Oo exp < — lUO z r > — U.loZW — U.l ( Oy lOr In ^ iO.D uev 

^\ Vl + 9.21(iy + 3.2G)/ J ^ y p- 
1.01 lor 'I], > GeV 




fen 


0.20463 tanh{ - 6.2370(2/ + 2.2)} - 0.16362(2/ + 1.6878)^ + . . . 






. . . + 3.5183 X 10-4(2/ + 9.6400)* 


a 


61 


r / 2/ + 2.0154 

1.6537 + 3.8530 exp <^ -3.2027 L ^ „ L^.,„, \ 

\ V1.0 + 0.62779(2/ + 2.0154)/ J 




fe2 


-10.722 - 0.082672 tanli{ - 1.8879(2/ + 2.1)} + 1.4895 x Vi^^{y + 256. 63)^ 


Diffrac 


fes 


r / ?/ + 1 9877 \ 2 1 

-0.023752- 0.51734 exp -3.3087 — — ^ 7tt7;==7 ? 

VI.O + 0.40300(2/+ 1.9877)/ J 




60, ... ,63 


for Tp < 5.52 GeV 




64 


0.94921 + 0.12280(2/ + 2.9)^ - 7.1585 x 10"*(2/ + 2.9)* + 0.52130 fogio {y + 2.9} 




65 


-4.2295 - 1.0025 tanh {9.0733(2/+ 1.9)} - 0.11452(2/ - 62.382) 




66 


1.4862 + 0.99544y - 0.0427632/^ - 0.00400652/^ + 0.00579872/* 




67 


6.2629 + 6.9517 tanh { - 0.36480(2/ + 2.1)} - 0.026033(y - 2.8542) 


o 




0.37790 exp j -56.826 f — — ^^'L 1 - 0.059458 + 0.00965832/^ 
^\ V1.O + 1.5221(2/ + 2.9537)/ j " 


o 


di 


-5.5135- 3.39881/ 


a; 


d2 


-7.1209- 7.1850 tanh{30.801(2/ + 2.1)} +0.35108y 




da 


-6.7841 - 4.83851/ - 0.915232/^ 




^4 


-134.03 - 139.631/ - 48.3162/^ - 5.5526^^ 



Table 2.5: Kamae et al. parameters describing electron spectra for arbitrary proton energy. 
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Parameterization for Positrons 




Parameters 


Formulae as functions of the proton kinetic energy {y = logjQ(rp)) in TeV. 




Oo 


-0.79606(2/ + 3.3) + 7.7496(i/ + 3.3)^ - 3.9326(j/ + 3.3)^ + . . . 




. . . + 0.80202(2/ + 3.3)* - 0.054994()7 + 3.3)'' 




ai 


6.7943 X 10"'^' - 3.5345 x IQ-^'y + 6.0927 x lO^^i/^ + . . . 




. . . + 2.0219 x lO^V + 5.1005 x 10"**/ - 4.2622 x lO^V 


iction 




44.827 81.3781ogio {0.027733(2/ + 3.5)1 _,_ g 44^7)2 




Q3 


0.52010 + 0.593363/ + 0.012032j/2 - 0.00642422/=^ 


D 


04 


2.1361(2/ + 3.32) + 1.8514(2/ + 3.32)^ - 0.47872(2/ + 3.32)^ + . . . 


o 

z: 


. . . + 0.0032043(2/ + 3.32)" + 0.0082955(t/ + 3.32)-' 

1.0845 X 10"'^ + 1.4336 x lO^^'logn, {0.0077255(y + 4.3)} + ... 




05 


1.3018 X 10-* s 

(, + 4.8188)^ +9-^^"^^^" 








"6 


-267.74 + 14.175 logio {0.35391('i/ + 3.4)} + 64.6692/ - 7.7036?/^ 




07 


1 9869 X 10* 

138.26 539.841ogio (0.12467(2/ + 3.9)} (^y 'j ^^'^/^^ 




as 


-0.14707+0.401352/ + 0.00398992/^ - 0.00166022/^ 




' KifJ 


f „„„/ V + 3.25 \2l 
2.22exp |-98.9(^ ^ 10.4(2/ + 3.25)) | ^ 

1.0 for Tp > 5.52 GeV 




bo 


29.192 tanh { - 0.37879(2/ + 2.2) } - 3.2196(y + 0.67500)= + 0.0036687(y + 9.0824)* 




bi 


142.97+ 147.86 expi 0.37194 f 2/ + 1-8781 \21 
^\ 11.0 + 3.8389(2/ + 1.8781)/ J 




62 


-14.487- 4.2223 tanh { - 13.546(2/ + 2.2)} + 1.6988 x 10-*(2/ + 234.65)= 


:actior 


63 


f / «; + 18194 \2~1 

-0.0036974 - 0.41976 exp-^ -6.1527(— — ^ I } 

^\ V1.O + 0.99946(2/+ 1.8194)7 / 


te 


60, ■ • • , 63 


for Tp < 11.05 GeV 


Q 


64 


1.8108 + 0.18545(2/ + 2.95)^ - 0.0020049(2/ + 2.9)* + . . . 




. . . + 0.85084 exp|-14.987(2/ + 2.29 - 0.18967(a; + 2.29))^| 






9 ndCld — n 51 54StaTiVi /9 9758('5/ -1-1 Ql\ — fl n'?5nnQ('5/ — fi f>555^ 




'-'6 


1 5258 + 1 01 32?/ — n 0643S8J/2 — 0040209?/^ — 0082772?/* 




67 


3.0551 + 3.5240tanh{ - 0.36739(2/ + 2.1)} - 0.13382(2/ - 2.7718)= 




Co 


f / u + 3.1272 9.6984 „ 

2.9841 exp<^ -67.857( „ ^ „ ^^„„, , _.Q^o^o^ K 6.5855 + 0.412562/= 

\1.0 + 0.22831(2/ + 3.1272)/ J 2/ 




Cl 


f; S97R 4- 99'^^^7^ ^ 1 zif^'?ni»2 


< 


C2 


0901 zL^si ^^anii / ^ ns'^nf^K 4- 9 1 ^\ 4- n At^'i^Oti 

— o.uzyi — u.4;0oi lami ^ u.uooui^2/ ^--^ / f ^ u.^uoozc/ 


*^3 


59300 + 36093?; 




C4 


0.77368 + 0.447762/ + 0.0564092/= 




do 


^■9^««^^p{-5'^-^^<i.0 + L2892(l + 2.9485)) }" 0-23720 + 0.0413152/= 


o 


di 


-4.9866-3.14352/ 


a; 


d2 


-7.0550 - 7.2165 tanh(31.033(2/ + 2.1)) + 0.385412/ 




d; 


-2.8915 - 2.14952/ - 0.45 06?/= 




di 


-1.2970 - 0.139472/ - 0.411972/= - 0.10641?/^ 



Table 2.6: Kamae et al. parameters describing positron spectra for arbitrary proton energy. 
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Due to the low statistics at very low energy, Badhwar's and Tan's parametrizations tend 
to produce unphysical distributions for proton kinetic energies below 6 GeV. Nevertheless, 
the total inclusive CS - the integrated version of those - are still in agreetment with the 
available experimental data. To fix this undesiderable feature, both parameterizations are 
patched by doing a smooth transition from 3 GeV until 7 GeV with the Stecker's model. 
Let's clarify that Kamae's parameterization also includes that feature, but considering 
more resonances. 

Moreover, for proton energies above 100 GeV, Badhwar's parameterization becomes 
unstable specially for the electron CS case. 



2.3 Production in annihilation processes 

Another phenomenon where positron and electron production takes place is in annihilation 
processes. Those particles can be produced directly and/or from the subsequent decays 
of other particles. For example, a muon eventualy produce an electron. As well, similar 
situation would happen when a quark hadronizes. 

An annihilation process is described by specific rules, which lay on some particle 
physics theory, for example, Supersymmetric models describe the Neutralino DM physics 
[551 EU]. Other types of annihilations are not out of that situation, in fact, the 
annihilation can be labeled in terms of outgoing particles (Figure I2.3|l . 



In that sense, the positron (electron) multiplicity distribution (MD), i.e. the number of 
positrons (electrons) per unit of energy produced in a single annihilation event, is: 

=5:BRt«^.) (^) . (2.42) 

and it depends on branching ratios of intermediate states i, 

BR(xx-^) = ^^^^^. (2.43) 

Ctotal 

Note that those are described by the theory which describes the annihilation process. 
Also, it is needed the MD that the z-state can generate, which is generally described by 
SM processes. 



In astrophysical scenarios, it is expected that DM particles move at non-relativistic 
speeds, i.e. ECM should be closely equal to the double of DM mass. Related to this, anni- 
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Lineros' PhD thesis [2008] 



LinEros' PhD thesis (2008) 



Badhwar et al. parameterization 




10 ^ 10" 
energy [GeV] 



Figure 2.2: Plots of positron and electron inclusive CS versus energy for proton energies of 2, 10 
and 35 GeV are shown. Each plot is obtained by using Badhwar et al. [48], Tan and 
Ng [50], and Kamae et al. [55] parameterizations. Also the relative difference respect 
to Tan and Ng solution shows a mean difference around 15% and 25% for Badhwar 
et al. and Kamae et al. cases respectively. 
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Linoms' PhD thesis (2008) 




Figure 2.3: Annihilation event which produces positrons and/or electrons is labeled by an in- 



hilation products should be composed by lighter particles which will correspond to SM ones. 

In the case of DM annihilations, intermediate states are described by the production 
of two SM particles which keep conserved quantities as in the annihilation; those are 
electrically neutral, colorless and unpolarized processes. 

The MD's can be calculated from analytical expressions in the simplest cases, but in 
complex ones, it is needed to use more sophisticated methods. 

2.3.1 Generation of Multiplicity Distributions 

Multiplicity Distributions for positrons and electrons are classified through their total 
electrical charge. As well, each one is function of ECM and kinetic energy of the target 
particle, in our case positrons and electrons. 

The method to generate them employed here is based mainly on the Lund's PYTHIA 
montecarlo generator [61j . which is used for simulating decay chains and hadronization 
processes. With PYTHIA it is possible to generate a basic set of intermediate states 
(Table I2.7|l . even though, extended set of states can be composed from the PYTHIA's 
ones. 

Generation from PYTHIA 

PYTHIA is used to simulate decay chains and hadronization processes from an initial 
(intermediate) state. Not all possible PYTHIA's initial states may produce positrons 
and electrons. Table [27f] lists all processes; however, some states as (e+e^), (/i+/i^) and 



termediate state i. 
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Intermediate state 


Charge 


Lcptons 


Quarks 


Gauge Bosons 


+ 1 


(z^ee+) (1/^^+) (VrT^) 


(ud) (cs) (t5) 


{-/W+) {ZW+) 





(e-e+) (r-r+) 


(uu) (dd) (cc) (ss) (&&) (ti) 


(gg) (ZZ) {W-W+) 


-1 


(e~Pe) (/^^l^Ai) (t^I^t) 


(du) (sc) (6f) 


(w-j) (w-z) 



Table 2.7: Intermediate states generated with PYTHIA for positron and electron MD's. 



its combinations with neutrinos will produce MD's which are obtained from analytical 
expressions (See Appendix |b|. 

The method for MD generation is based on histogram creation from simulated events. 
Histograms correspond to number of positrons (electrons) per kinetic energy bin, where 
each bin follows a Poisson distribution. Moreover, the number of total positrons (electrons) 
produced is A^part — 10^ to reduce as maximum as possible statistical uncertainties. As 
well, to reduce even more possible uncertainties each histogram is done using the Method 
of Average (See Appendix |e|) 

An extra feature, which is not included in the original PYTHIA program, is the 
addition of muon's polarization effect (See Appendix |b|) . This is performed by a selection of 
muons - in generated events - which have been produced from spin-0 particles, as mesons 
and kaons (See Appendix |c] and Appendix |d]). This effect produces differences up to 10% 
between positron and electron MD's at the same kinetic energy. 

Each state has been simulated for different ECM. The energy range covered is from the 
threshold energy till 20 TeV. In that way, most of possible interesting DM mass range is 
covered. Although there are special values of ECM, 

Vs"spcciai = "^T, 2mr, rub, 2mb, mw, mz, 2mw, rut, 2mz, 2mt, (2.44) 
which are important keys in composed states production. 



We generated 60 MD's - in the range previously described - for each intermediate state. 
Furthermore, interpolations are performed for the positron (electron) energy and ECM. 
Interpolation variables are x and t: 

variable x usually goes from to i. The interpolation is performed in two steps, where the 
first one is on x and then on r. Let's specify that instead of interpolating the MD values 
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Energy [GeV] Energy [GeV] 

Figure 2.4: Multiplicity distribution for positrons and electrons obtained from PYTHIA. Plots 
made for DM masses of 100 GeV and 1 TeV for W^W' , t'^t~ and bb intermediate 
states (channels) and the tt channels with masses of 200 GeV and 1 TeV. 



we interpolate the logarithm of this. In that way, the interpolation is stable and correctly 
reproduces MD's for intermediate values of x and r. 



Standard Model composed states 

Montecarlo generator programs, as PYTHIA, are useful to construct MD's, however, there 
are still states where those are not an efficient choice. For those solution is to 

compound them by using already known MD's. 

In the SM particles context, MD related to higgs particles (Table I2.8|l should be 
composed because their mass value is still unknown. The composition is based on higgs 2 
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SM intermediate state 


Charge 


Higgs-Higgs 


Higgs-Gauge 


+ 1 




{hW+) 





(hh) 


{hi) (hZ) 


-1 




(W-h) 



Table 2.8: Composed states for SM. 



body decay modes, which are dominant. 



Decay widths, in which higgs goes to a couple of leptons or quarks, are computed from 
SM Feynman rules: 



nh - //) - iT^^^^h 1 - ^ , (2.46) 
where Nc is the number of color and it takes values 1 for leptons and 3 for quarks. 



Higgs can also decay into gauge bosons: 



T{h^ZZ) = ^^yr^^fl-xz + Jxl) (2.47) 
1287rm^ V 4 y 

Tih^W+W-) = ^I!!^^T^^(l-:cw + ^xl) (2.48) 



where xz/w = 4m|/iy/m^. 



There are extra decays modes, such as /i — *■ 77, jZ and gg. However those are one-loop 
processes, and can be safely neglected in most of cases [52] . 

The positron (electron) MD for a single higgs at rest is composed by: 



{ee,Vs=mh), (2.49) 



where i-states correspond to II, qq, ZZ or any allowed decay mode. To compose the higgs' 
MD the ECM should be equal to the higgs mass. 

Nevertheless, it is almost improbable to produce a higgs at rest. To generalize the situ- 
ation, the higgs MD should be boosted to a reference frame where its energy matches the 
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production energy. A proper procedure is to use the boosted spectra formahsm (Appendix [a|), 
where for positrons and electrons is: 

where 7 = Sh/mh and the integration hmits are: 

t = min(emax,(7 + /3)£e)), (2.51) 
r = (7-/3)ee. 

Note that ^max is the maximum allowable energy in the higgs rest frame and it also helps 
to calculate the maximum energy in the boosted frame, 

Ee.max = (7 + /3)Cmax. (2.52) 



Two-particle states MD are created as direct sum of former single-particle MD: 

each single-particle state is set to conserve energy and momentum: 

SA^ TTr^ — - and eb = tTt^ -■ 2.54 



For example, the state hh is the simplest to be computed, as it is formed by two identical 
particles. Then, the positron (electron) MD is: 

|i)^__,e.,Vr) = 2x(Ji)__,..,.,. = f), (2.55) 

that is just two times a single higgs MD. 

A different situation happens for states hZ and hW^. Those states are formed by 2 
kind of particles. The Z and single states are obtained from two-particle states jZ 
and jW^ by taking advance of that photons do not produce any positron or electron. The 
single-particle state and the two-particle one are related as follows: 

(j-j i£e,£z/w) ^ [j^] {ecVs' = ez/w +Pz/w), (2-56) 

where -fZ and jW^ are computed with PYTHIA. 
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Lineros' PliD thesis C2008) Lineros' PhD thesis (2008) 




10 ^ 10° lo' 10^ 10^ 10 ^ 10° lo' 10^ 10^ 

Energy [GeV] Energy [GeV] 



Figure 2.5: Multiplicity distributions versus energy for SM higgs-higgs and higgs-Z composed 
states. Examples with DM masses of 200 GeV and 1 TeV are shown. 



The final step is to combine them by direct addition: 

(^) (e,,V^)=(^) (se,eh)+(^) iee,ez/w), (2.57) 



where 



^ 2^ '^^^^ ^ 2^ ' ^ ^ ^ 

Let's emphasize that this method works well for most of the states which do not need 
quarks as single-particle MD's. 



For example, Figure 12 . 51 shows MD's for composed states hh and hZ, assuming a 
higgs mass value of 125 GeV. The hh state produces more positrons and electron than 
hZ state. That is related to how they couple to fermions, especially quarks, which are 
major contributors to positron and electron production. Higgs bosons strength interaction 
is proportional to fermion mass, instead of Z bosons that couple in the same way with 
all fermions. Then, it is expected that higgs bosons decay mainly in bb pairs, where 
hadronization processes have a strong effect in the production. 



Two Higgs Doublet Model composed states 

As an alternative for the SM higgs sector is the Two Higgs Doublet Model (TDHM) which 
includes SM fields and two higgs isospin doublet - instead of one. Those may have equal or 
opposite isospin charge. An advantage of this model is to describe a richer higgs sector as 
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supersymmetric models do. 



The observable particles are increased in number, where it appears three neutral and one 
charged higgs fields. This makes possible new MD's (Table [2~9|) . Following similar procedure 
as SM cases, those are generated in terms of two body decay modes described by the model. 

Two body decay widths are easily calculated from model's Feynman rules 62] . Charged 
higgs decay are mainly dominated by decays into heavy fermions: 

327rm^m'^+ 

(m^+ — — ni^){mfj tan P + cot /3) — Amf^m^ 1 , 

where Nc is the number of colors, tan/? is the ratio between vacuum expectation values of 
two doublets and A^/^ is a kinematical factor, 

A"^/^(mo, TOi, 7712) = \J {ni\ + m| — m|) — Arnlm^ . (2.60) 

Also there is another decay mode into a and a lightest neutral higgs: 

647r7re^+ 



2 , 2 \ , 



where a is the mixing angle among neutral higgs fields. Furthermore, this mode is less 
important than previous ones because some regions in the parameter space produce 
physical configurations, for example charged higgs mass has to be bigger than the sum of 
lightest higgs and bottom quark masses. 

Decay widths for neutral higgs are quite similar to SM ones. These can be written in a 
compact way 62 : 

Ncg'^m'^dlm.h ( ^ Amj"^^ 



r(/.^»)=,yT'°^?Jl-^)'. (2.63) 



327r777^ cos^ 13 
where dh and Ch are defined as: 

- sin a h ^ Hq cos a h — Hq 

dh = { cos a h = ho , Ch = i sin a h = ho , (2-64) 

cos f3 h = Ao I — sin (3 h — Ao 
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TDHM intermediate state 


Charge 


Higgs-Higgs 


Higgs-Gauge 


+ 1 


{hoH+) {HoH+) {AoH+) 


/ 7 11' _A_ \ / -rr II' _|_ \ / A II' _J_ \ 

{hQW+) {HoW+) {AaW+) 
(7i/+) iZH+) 





(hoho) (HoHo) (AoAo) 
(hoAo) (HoAo) (H+H-) 


{W+H-) {H+W-) (hoj) (i/o7) 
(Ao7) (hoZ) (HoZ) (AoZ) 


-1 


(hoH-) (HoH-) (AoH-) 


(hoW-) (HoW-) (AoW-) 
(7i7-) (ZH-) 



Table 2.9: Composed states for TDHM. 



and the power index p is: 

_ f 3/2 ho,L 
1/2 h = Ao 



With most important decay widths aheady calculated, branching ratios are also known. 
Those help to weight each MD for composing a single higgs MD: 

where h could be ho, Hq, Aq or H^. Note that charged intermediate states are specially 
useful for computing charged higgs MD, which increase accuracy of the final MD. 

Depending on the DM candidate, different kind of higgs particles are involved. Compar- 
ing SM with neutral TDHM higgs fields at equivalent production conditions (Figure r2.6|) . 
we see how SM higgs produces more positrons and electrons, instead ho and Aq higgses, 
for higher values of tan /3. This is related to neutral higgs field couplings, which depends 
directly of parameters (3 and a. 



31 



CHAPTER 2. POSITRON AND ELECTRON PRODUCTION 



Lineros' PhD thesis C20Q8) Lineros' PhD thesis (2008) 




Energy [GeV] Energy [GeV] 

Figure 2.6: Multiplicity distribution for positrons and electrons versus energy. Comparison be- 
tween SM higgs and TDHM neutral higgs with equal mass value for a DM mass of 
1 TeV. 
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Positron and Electron 
propagation in the galaxy 

The propagation of positrons and electrons, and in general of any CR, in 
the Galaxy can be a hard problem to solve. There are many possible ways 
to deal with the propagation. The Two-Zone Propagation Model provides a 
good physical approach to describe - in a nice way - the travel between the 
CR source and the Solar System borders. Up to this point, the influence of 
the Sun becomes strong enough that CR need to be specially treated. 

3.1 Overview 

In the travel across the Galaxy, cosmic rays are affected by many processes. In the low- 
energy range, one of the most important process is magnetic diffusion, which is produced 
by random magnetic regions that fill all the Galactic surrounding. These magnetic regions 
affect the cosmic ray propagation in such a way that makes impossible to trace back a 
cosmic ray up to its source. As well, cosmic rays interact with the gas and other particles 
present in the interstellar space. This could make cosmic rays loose (or gain) energy or to 
produce other types of cosmic rays. 

Most of the sources of Galactic cosmic rays are in the Galactic Plane (CP). However, 
sources can be also located outside of it. In this way, a model to describe CR propagation 
should be able to represent all this possibilities. 

Apart of the Galactic-scale effects, which rule the propagation, we cannot forget the 
effect associated to Solar activity. The Sun has 11-years cycles during which it increase 
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and decrease periodically its activity. One of its manifestation is the increment of the Solar 
Wind flux, which acts on CR by pushing them away from the Solar System and reducing 
their energy. 

In this chapter, the Two-Zone Propagation Model is described and solved for the case 
of positrons and electrons. General solutions of the transport equation are explained. As 
well, the model's space of parameters is discussed and the solar modulation problem and 
standard method to model it are also reviewed. 



3.2 Two— Zone Propagation Model 

This is a model for studying CR propagation in the Milky Way In general terms, 

it is composed by a Propagation Zone (PZ) which demarcates the region where CR 
propagates and by the Transport Equation (TE), which models the physics of propagation. 

The PZ is composed by two cylinders centered at the Galactic Center (Figure 13. 
Both cylinder has a common radius equal to the galactic one (R = 20 kpc). However, their 
thickness are rather different. 

The thick cylinder has a height of 2L and fills all the PZ. Its height L is related to 
how much magnetic fields extend. CR measurement have constrainted its values to a range 
that goes from 1 to 20 kpc [M]. Moreover, in the thick cylinder, processes as magnetic 
diffusion and energy losses - related to interaction with magnetic field and light - take place. 

The second cylinder is a thin disk with height equal to 2hz, where hz — 100 pc. 
The interstellar medium, cosmic-ray sources and interactions, as energy losses related to 
ionization or Bremsstrahlung, are contained inside the thin disk. Also, CR reacceleration 
processes related to supernova explosions and Shockwaves take place there. 

In the PZ, the Solar System is placed inside the thin disk at 8.5 kpc from the GC and 
lays in the Galactic Plane. 



3.3 The Transport Equation 

The Transport Equation describes the physical effects involved in the CR propagation inside 
the PZ. In general terms, TE is based on a continuity equation for the CR density per unit 
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Lineros' PhD thesis (2008) 



R 



L 




Figure 3.1: Propagation Zone geometry for the Milky Way. A cylinder of radius R {— 20 kpc) 
with thickness of 2 L delimits the region where CR propagate. A small cylinder with 
same radius but with thickness 2 hz (= 200 pc) models the Galactic Plane. The Solar 
System is placed at the Galactic Plane with distance tq = 8.5 kpc from the Galactic 
Center. 

of energy (?/') that written in terms of currents is simply: 

dj'{t,x,e;^,dilj) = s{t,x,e) , (3.1) 
where the currents J'{t,x,e;Tp,d'4') denotes CR displacement in time, space and energy, 

dj^ + dj, +dJt^ s{t,ii,e) , (3.2) 
where s(i,x, e) are sources (or sink) which are independent of cosmic-ray evolution. 

Diffusion due to random magnetic field and interaction with Galactic Wind (GW) are 
contained into the spatial current: 

Jx = -Dad^i; + Vi/- , (3.3) 

where Dq is the diffusion coefficient and V is the vector field of the GW. Usually diffusion is 
treated as a homogeneous function in space, dependent on the CR rigidity and the Lorentz 
(3 factor [Ml 122] as, 

Do{e) =KoP X (^^^ , (3.4) 
where TZq is a rigidity scale, usually 1 GV. fn the case of ultrarelativistic particles it is just: 

Doie)=Koe^, (3.5) 

where e = e/eo with Eq — I GeV. 
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There are some models to explain the form of the Galactic Wind field [Ml EZl [Ml [M] > 
however it is still not well determined. The choice is to assume a GW field perpendicular 
to the Galactic Plane, 

V(x) = V,{z)z . (3.6) 

Furthermore, it is assumed a constant behavior when we are out of the Galactic Plane, 

Vc{z) = K sign(z) . (3.7) 

Nevertheless, this choice may be not unique. In the work of Bloemen et al. f7(T, they have 
used a lineal profile, Vc{z) — 3 Vq z, for modeling the GW. This profile has the advantage 
that in some cases, the TE can be solved analytically. 



The energy-current (J7e) is related to energy losses and gains. 



de 

= -0 = Vss/gain(e) • (3.8) 



In the galactic environment, CR are affected by processes like bremsstrahlung and 
ionization of IS gas, synchrotron radiation and inverse Compton scattering with photons 
and magnetic fields. As well, GW affects them producing adiabatic losses, i.e. CR are 
cooled due to the expansion that GW drift produces in them, but just in zones where V 
variates spatially. On the other hand, CR are also allowed to gain energy, that is possible 
through reacceleration processes. 

The energy losses due to inverse Compton scattering and synchrotron radiation are 
calculated from the total radiation power 71] : 

where is the photon density. A typical value for Wph is 1 eV/ cm"^, which is a nominal value 
to include starlight, infrared and microwave radiation J2. . In the case of ultrarelativistic 
positrons and electrons, the energy loss term is: 

6synch(e) , (3.10) 

te £o 

where Eq = I GeV and te = 10^^ sec is a time scale for the synchrotron energy loss. 

Another process to consider is the adiabatic cooling due to the GW interaction. The 
energy loss term mainly depends of the divergence of the GW vector field, 

badiab(£,x)^- ^ - , (3.11) 

6 £ 
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that is reduced, in the TZPM context, to: 
&adiab(e,x) = 



i dzVc{z) e inside GP 

(3.12) 
outside GP 



for ultrarelativistic positrons and electron. Note that adiabatic cooHng becomes important 
inside the GP: this puts in evidence how the inner structure of GW may take an important 
role in the propagation. 

Other processes like bremsstrahlung, ionization and Coulomb losses, have contributions 
in the range above ~1 GeV smaller than the inverse Compton scattering and synchrotron 
radiation |73| . For example, the ionization energy losses for ultrarelativistic positrons (elec- 
trons) in neutral hydrogen and helium are: 

^ion(e) = -A'ion ^ Z^n^ (31og(7) + Ci) (3.13) 

i=H,Hc 

where Zi is the atomic number, Ui is the gas number density and the constant 
Kion = 7.484 X 10^^^ GeV cm'^/s. The parameters Ci a-re related to the ionization potential, 
where the numerical values are: Ch = 20.46 and Cho = 19.27. Notice that for a electron 
energy of 10 GeV. We got that |6ion| ^ 10^^® GeV/sec, which is two order of magnitud 
smaller than the same case but for energy losses due to synchroton radiation and inverse 
compton scattering |&synch| ^ 10^^^ GeV/sec. The same is illustrated in Figure 15. 21 where 
time scale for those are shown. Smaller time scale interaction means bigger energy losses. 

On the other hand, CR may gain energy through interactions with Shockwaves. In the 
TE, this is included as an energy-gain term and an energy-diffusion term: 

Js.rcac = i^A'ee V' " /S'i^ee 0,^ (3.14) 
£ 

where K^e is the energy-diffusion coefficient, 

4 y„2p2 

= 3S{A-S^)iA-S) D^) ' ^^-^^^ 
which depends on the Alfven velocity (Va) related to velocity of disturbances in the hy- 
drodynamical plasma present in the ISM [75] . These expressions are simpler in the case of 
ultrarelativistic particles: 

Je.rc.c = ^K,,i;-K,,d,^ and i^,, ^ }^ ^2-5 _ (3^^g) 

Note that the inclusion of reacceleration processes in the TE produces a big change into 
the nature of the differential equation and also in the strategies to solve it [SS] . 
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Figure 3.2: Time scale related to energy loss processes for positrons and electrons as reported 
in [73]. In the energy range above 1 GeV (orange line), energy losses are dominated 
by inverse Compton scattering and synchrotron radiation processes |74) . 



Time evolution of CR are important just in the case of a transient source. But for typical 
cases, like those discussed in the next chapter, time evolution can be safely neglected. In 
the energy range above 1 GeV some of energy related processes can be safely neglected, as 
well. 

In this context, the positron and electron transport equation (PETE) is reduced to: 

-Koe'V^^-^(—A ^s. (3.17) 



de \te 

It includes the main processes (magnetic diffusion and inverse Compton scattering and 
synchrotron losses) that affect positrons and electron during their propagation inside the PZ. 

The cosmic-ray flux is the quantity that experiments are able to measure. The solution 
of PETE gives just the information about the cosmic-ray density. The cosmic-ray flux, up 
to the Solar System borders, depends directly on the cosmic-ray density, 

^(e) = g^(x0,£). (3.18) 

Let us warn that the evolution from the Solar System borders until the Earth can not be 
done through solving PETE, because the influence of Sun's activity makes this approach 
not reliable. Instead of this, the problem of propagation inside the Solar System can be 
treated in a simple way by modulating the flux, as will be discussed later in the chapter. 
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3.4 PETE solutions 

A general solution of PETE is found using the Green function method. The Green function 
(G) is the solution for a point-like source in space and energy, 

- Ale'V^G - ^ (e^G) = te S^^ - x,)(5(e - e,) , (3.19) 
where Aq = KoTg and formally G = G(x, Xg, e, Cg). 

Furthermore, this type of differential equations can be solved by the method of sepa- 
ration of variables. The Green function and the Dirac delta on space can be expressed in 
terms of solutions of the Helmholtz equation (Xg), 

G(x,x„e,e,) = ^03(e)xg(x), (3.20) 
g 

^3(x-x,) = ^xJ(x«)Xs(x), (3.21) 
g 

where Xg satisfies: 

+ g\g{^) = . (3.22) 

After the substitution into PETE, we found that the equation is just a set of ordinary 
differential equations labeled by the eigenvalue g: 

Ale'g^^e) - ^ {e'M^)) = xl{^s)5{^ ~ e.) • (3.23) 
When e 7^ eg, all the equations become homogeneous and with solution of the form: 

<^,(e) = 5exp(.g2Ag|^) . (3.24) 

A restriction arises when the non-homogeneous equations are integrated around e^, 

.^+(e,) -</-;(£,) = -^xJ(x.), (3.25) 

which links the solutions for energies above and below of e^. Furthermore, we expect that 
particles are not allowed to gain energy because our system describes propagation with 
energy losses only. This implies that: 

(/.+ (£.) = =^ ^-(e,)^I|4(x,), (3.26) 

which fixes the initial conditions for all solutions, we obtain in that way that the solution 
of non-homogeneous equations is: 

U^) = ^Xj(x«) exp (-1 , (3.27) 
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where 



Ad = 2AoW ' (3.28) 

is called the diffusion length. And finally, the Green function is: 

G(x, x„ e, e,) = -| G(x, x„ Ad) , (3.29) 
which is proportional to the tilded Green function: 

G(x,x«,AD) = 5^xJ(x«)Xs(x)exp(-i f/A^,") . (3.30) 



4 

Notice it depends on energies through Ad- This function can also be described as: 

G(x,x„Ad) = (x«|w(Ad)|x) , (3.31) 

where w{Xu) = exp(— | g^^u) is the weight for each state Xg- As well, this function 
satisfies: 

lim G(x,x„Ad) = <53(x-x«) , (3.32) 

Ad— »0 

which describes a zero-distance propagation limit. In other words, particles observed with 
e = £s at X cannot come from anywhere, except from the source at same position. 

The general solution of PETE is obtained by convolution between the source term and 
the Green function: 

ip{x, = ^ j d^Xg q{xs, eg) G(x, x^, e, e^) . (3.33) 



3.4.1 Solution in free space 

Solutions of Helmholtz equation in a three dimensional free space are: 

^'^ " (27r)3/2 ^""P k • x) , (3.34) 

where the eigenvalue k is the composition of the eigenvalues that correspond to each spatial 
dimension, 

k ^— kj. ~\~ ky ~\~ k^ . (3.35) 
In this case, the tilded Green function is obtained from: 

Gfree(x,Xs, Ad) = J d^k exp ^ik • (x - x^) - ^/c^A^j , (3.36) 
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note that this is separable into three identical gaussian integrals of the form: 
J dfc exp (-fc^a^ + i/c6) = exp ^-^^^ , 



(3.37) 



one for each dimension. Each of them produces its own tilded Green functions, which are 
independent from the others, 

- exp 



Gfreei^, Xs,\d) ^ —^=rr- exp "2 . (3.38) 



This independence helps to compute the tilded Green function in three dimensions: 

Gf,t(x,x„AD) = Gl^,,ix,Xs,XD)GlLiy,ysAB)GlLiz,Zs,XB) (3.39) 

1 / (x-x^\ 

At this point, it becomes clearer the name given to Ad because it takes the place of the 
diffusion length in the standard diffusion theory. Note that this procedure produces the 
same results as reported by Baltz and Edsjo [ID]. 

Let us denote that Ad is a key quantity in order to understand the propagation. In 
Figure 13.31 we see how the source distance modifies the value of Gj^f^^ , but on the other 
hand the maximum value is reached when Ad is of the order of the distance between the 
observer and the source. 



As seen before, the solution of PETE depends on the convolution between the green 
function and the source term (Equation 13. 33|l . By inspecting this, we note that the leading 
contributions to the solution come from a sphere of radius Ad centered at x. 



3.4.2 Solution with boundary condition 

The TZPM considers propagation to occur in a finite volume, where CR may escape if 
they arrive close to the boundaries. 

A first case is to consider boundary conditions on the vertical axis, 

Tjj{x,y,z = ±L) = . (3.40) 

As before, there is a orthonormal basis for x, y dimensions, which is described by a con- 
tinuous 2D Fourier space - Equation 13.341 However, to satisfy the boundary condition on 
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Lineros' PhD thesis (2008) 



Zg = 0.1 kpc 
Zj. = 1 kpc 
Zj. = 10 kpc 




Figure 3.3: Gfj.cc versus \t>I Zs. The green function are smaller for farther sources. The maximum 
value is reached for Ad of the order of the distance to the source. 



z, a discrete version is required: 



where the eigenvalues fcz are: 



sa\(kz.n z) n : even 
cos(fcz,n z) n : odd 



nn 

~ n E £i 



(3.41) 



(3.42) 



In the tilded Green function, two parts are identified. The first is similar to the free case 
and the second part satisfies the boundary conditions: 



G(x,Xs,Ad) = |Gfree(a;,a;s, Ad) Gfree(y,2/s, Ad)| X I](z,0s,Ad), 
where S is calculated as a superposition of Fourier modes: 

^f{z,Zs,Xu) = ^ X! sin(fc^_„0s)sin(fc2,„2;)exp f --fcf „Ad j 

n:cvcn ^ ' 

1 °° /I 

+ — ^ cos(A:^,„Z5)cos(fc^,„z)exp ( --fcf „Ad 



odd 



(3.43) 



(3.44) 



A different way to compute it is using the method of Image Charges (IC), where S is 
the superposition of many free-space Green functions [40] : 



Sic(z,Zs, Ad) 



n—oo 

E 



(-l)"Gfree(^,^.,n,AD) 



(3.45) 
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10 ^ 10 ^ 10° 10^ 10^ 



Figure 3.4: E function versus Xu/L. The IC and Fourier method E are calculated with few 
terms to contrast those with an exact solution. It is shown that the IC method (red 
line) converges faster for Ad values smaller than L instead that the Fourier one (gree 
line) which converges for Ad bigger that L. 

where z^.n denotes image charges positions, 

^,,„ = 2in+(-l)"z, , (3.46) 
needed to satisfy boundary conditions. 

Both methods are completely equivalent from the theoretical point of view. However, 
depending on the value of Ad/^, the rates of convergence are different (Figure 13 . 4|) . We 
found two regimes in which each method becomes the most efficient option: 

i) When Ad < L, the IC method is more efficient than the Fourier one. The sum 
converges after few terms because the diffusion length (Ad) is small enough that 
the propagation cannot reach the boundaries. On the contrary, the Fourier method 
becomes highly inefficient because the weight function (Equation 13 .311) becomes unity 
for many Fourier modes. That affects directly the rate of convergence because we 
need to sum too many terms. 

ii) When Ad < L, the Fourier method is the best option because Fourier modes with 
low eigenvalue mainly contribute to the sum and we need just to consider few of them. 

In Figure 13.51 the effects of boundary conditions are shown. Boundary conditions 
induce the Green function E to decrease faster than the free-space one, when the 
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Figure 3.5: E and Gjf^^ functions versus \d/L. The vertical boundary conditions produces a 
green function which decreases faster than the free-space case when Ad grows. For 
diffusion distance smaller than L, both Green functions have same behavior. 



diffusion length is bigger that L. This is related to the leaking of particles which 
arrive at the boundaries. For the same reason, both Green functions are similar when dif- 
fusion lengths are small, particles do not have time to propagates and reach the boundaries. 



3.4.3 Axial symmetric solution 

The PETE can be also solved for special geometries. Galaxies, as the Milky Way, tends 
to have axial symmetry with respect to the galatic center and oriented along the angular 
momentum direction. In such scenarios, it is expected that CR sources are located 
according to the symmetry. 

In cylindrical coordinates an axial symmetric point-like source can be described as: 

q{r, z, e) - fclil S{z - z,) 5{e - e,) , (3.47) 
r 

and in this case the source corresponds to a ring. 

The function basis for the radial coordinate is based on Bessel functions, 

XkU^) = -fT^^.M^r.^^) ^ (3-48) 

where ao.i is the i-th zero of Bessel function Jq. This basis also includes boundary con- 
ditions, where Green function vanishes at a characteristic radius R. This condition implies 
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discrete eigenstates with eigenvalues: 



R ■ ^^-^^^ 



Following the general procedure, the tilded Green function for the radial coordinate is: 

T(r,r„ Ad) =2^2 -J— exp --k^^iX^ ■ (3.50) 

i=i [RJi{ao,i)) \ * / 

In the case with no radial boundary condition, this expression changes into a Bessel 
transform, which is not analj^ically solvable like in the cartesian case. 

The Grccm function with boundary conditions on the radial and vertical axis is just the 
product of the previously seen Green functions: 

Gcyi(r, rs, z, Zs, Ad) = T(r, r^, Ad) f^{z, Zs, Ad) • (3.51) 

Let us emphasize that T behaves as S when Ad is bigger than R. In a sense, particle 
leaking starts to be a dominant effect, producing a fast decrease in the Green function 
value. Generally L < R — vq, that means leaking in vertical axis will be even more 
dominant, because particles are closer to vertical boundaries than radial boundaries from 
Solar System position. 

In this case, the solution of PETE is: 

tp{r,z,e) = J des j r^drs J dzs q{rs,Zs,es) Gcyi{r,rs,z,Zs,XTi) , (3.52) 
where the region of integration is restricted to the zone inside the boundaries. 



3.5 The Halo Function 

The Halo Function (HF) is an dimensionless function which encodes the information about 
spatial dependence of the source term. It can be calculated when energy and spatial depen- 
dence in the source terms can be separated as: 

9(x,e) = KX/(x)x5(e). (3.53) 

Furthermore, /(x) is an dimensionless function which is normalized to the Solar System 
sources density, 

/(X0) = 1 , (3.54) 
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where the Solar System position is located at 

tq — 8.5 kpc and zq = kpc , (3.55) 

with respect to the Galactic Center. Let us remark that the normalization is just valid in 
cases of regular and non-vanishing distributions. 

Similar to PETE's general solution, the HF is calculated as follows: 

/(Ad) = y^^x, /(x,) G(xq,x„Ad) , (3.56) 

which is a convolution between the spatial part of source and the Green function. Also, 
HF naturally satisfies: 

lim 7(Ad) = 1 , (3.57) 
that is a direct consequence of the normalization of /(x) and the Equation 13. 321 

The main aim of the halo function is to be used for calculating the PETE's solution: 
^Q{e) = K^ de, g(e,) /(Ad) , (3.58) 



e 

where Ad has an implicit dependence on e and fEquation l3.28|l . 

In any case, HF are functions of one variable. This is an advantage that improves the 
calculation speed for several different cases. For example, many sources may share same 
spatial distribution but different spectral shapes, in this case we need to calculate only 
once the HF and just perform many convolutions in energy without recomputing the HF. 

A special case happens when source terms inject monoenergetic particles because the 
PETE's solutions become proportional to the HF, 

^Q{e) = K^I{Xu) . (3.59) 

But in general rule, any PETE's solutions depend directly on the spatial distribution of 
sources. 

An alternative approach to compute a HF instead of Equation 13 .561 is: 

7(Ad) = E 4 Xai^Q) (-\9^^d) , (3.60) 

which is obtained when the tilded Green function is decomposed in terms of solutions of 
the Helmholtz equation (Equation 13 .301) . The coefficients Ug, 



ag= I d-^x^ /(xs) Xgi^s) , (3.61) 
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are the projections of the spatial source distribution on the eigenstates of the Helmholtz 
equation. 

In analogy to the tilded Green functions computed with Fourier and IC methods, 
the two previous versions of HF have different convergence rates. The HF in terms of 
eigenstates (Equation I3.60|l converges faster for big values of Ad- On the contrary, HF - 
obtained by convolution - converge faster for small values of Ad, but just in the case when 
the tilded Green function corresponds to the one obtained with the IC method. 



3.5.1 Halo function examples 

A very simple example occurs for homogeneous sources distributed in all the space. In this 
case, we get: 

ifrcc(AD) = / d'x, G3d^(xQ,x,, Ad) = 1 , (3.62) 

J all space 

when the source is confined into a box centered on xq with edge lengths 2 l^, 2 ly and 2 1^ 
for each dimension, the correspondent HF is: 

7box(AD) - Erf (^^^ X Erf (^^^ x Erf (^-^^ , (3.63) 

where Erf(a;) is the Gauss error function. In this case, the disappearance of the source, 
beyond the box limits, starts to manifest, decreasing HF value, from diffusion lengths 
bigger than half of the minimum side of the box. Also note that when the box sides become 
bigger than Ad, ?box goes to /free- 

In the axial symmetric situation, we can calculate the HF for a cylinder with radius 
i?src (< R) and thickness of 2ft,2 (< 2L). The homogeneous source distribution makes easy 
to split the HF into an exclusively-radial and -vertical HF. This is possible when the source 
term is composed by functions that depend only of one spatial dimension. 

The radial part is described in the same way as a generic HF f Equation I3.60p . but in 
this case, it depends on Bessel functions: 

where /c^.z are the eigenvalues of the correspondent eigenstate f Equation 13. 49|) . In any case, 
the coefficients in the sum come from the radial integration over all the PZ, 



RJl{ciQ^i) Jq ' kr^i R Ji{kr,i R) 
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The vertical part is obtained for 2 cases. The first one is related to the Fourier method, 
which is similar to the previous examples: 

/vcrtical,F(AD) = J 2^ ^^^[^7 I • (3-66) 



A different approach is used to calculate the vertical part. In this case, the tilded Green 
function obtained with the IC method (Equation I3.45|l is used to perform a convolution 
with the spatial distribution, obtaining: 

/ve.ticaUc(AD) = ^ ^ Erf j - Erf . (3.67) 

As seen previously in the comparison between IC and Fourier methods, these two vertical 
parts present the same feature respect to convergence speed. The one based on the Fourier 
method converges faster for bigger values of Ad instead of /vertical, IC: and viceversa. 
Finally, the HF for a homogeneous source is: 

/axial(AD) — /radial(AD) X /vertical(AD) , (3.68) 

which depends on the radial and vertical parts that were already calculated. 

When the radial and vertical HF are combined to form /axial, the effect of vertical 
boundary conditions manifestates earlier ii L < R (Figure |3.7|) . This produce that particles 
can reach first the vertical bondaries than the radial one. Typically, R is fixed to the 
Galactic radius (20 kpc). This condition suggests that for Ad of the order of ^ 12 kpc the 
radial boundaries starts to affect the behaviour of HF. Nevertheless, this limit depends on 
the source ditribution. If the sources are close to GC (the farthest point to any boundary), 
their contribution to the HF are less modified than HF for sources close to the boundaries. 

Some HF, from sources distributed in such a way that can not be separated as the 
previous example, are calculated following the general prescription (Equation I3.60|) . In the 
case of axial symmetric sources the HF is: 



oc oc 



/axial(AD) = XI X! Mkr,t Tq) CXp ( -^.g,>^D ) - (3-69) 



i—1 n:odd 

where the coefficients are: 



R 



^''^ ^ L R"^ Jl{k R) J ^'^^ ' f{^^^') Jo{kr,i r) cos{ky^^ri z) . (3.70) 
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Figure 3.6: Plot of the radial (right) and vertical (left) parts of axial symmetric HF versus Ad- 
The influence of boundary conditions produces that each part - radial and vertical 
- vanishes for values of Ad larger than the closest boundary. For the radial part, the 
boundary effects start to dominate when Ad > 12 kpc. On the other hand, the same 
happends to the vertical part when Ad > 4 kpc. 



T. DelBhaye, R. Linpros, N. FornengQ. F. Doneto Sr P. Salali (2007) 



1,2 - with radial boundaries (Bessel) 

■ witiiout radial boundaries (Green) 




Diffusion length A„ = (4K(,t) '/^ [kpcl 

Figure 3.7: Halo function versus diffusion length for a homogeneous sources for values of L(= hz) 
of 3, 10 and 20 kpc {R = 20 kpc). The effect of boundary condition in the PZ is 
reflected as a fast decreasing of HF at big values of Ad- 
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label 


S 


Ko [kpcVMyr] 


L [kpc] 


Vc [km/ sec] 


Va [km/ sec] 


^B/C 


max 


0.46 


0.0765 


15 


5 


117.6 


39.98 


med 


0.70 


0.0112 


4 


12 


52.9 


25.68 


min 


0.85 


0.0016 


1 


13.5 


22.4 


39.02 



Table 3.1: Astrophysical parameters giving the maximal, median and minimal antiproton flux 
compatible with B/C analysis |77) . 



The case of small Ad is to calculate the HF neglecting radial boundary conditions: 

4xiai(AD) = / Tsdrs I dzs f{rs,Zs) J^ic{zq,Zs,Xd) ThccirQ,rs,XD) (3-71) 

which depends on the function Sic (because it is fast to converge) and Tfi-co, which is a 
tilded Green function similar to T (Equation l3.50l) but without radial boundary. This Green 
function has an analytical form, 



Tfree(r, r., Ad) = ^ exp ( -'-^-j^ Wo ) (3-72) 



where 

/^(x) = e-" Io{x) , (3.73) 

which depends of the modified Bessel function /o(a;) [75]. In the regime of big x values, 
that is equivalent to small Ad values, this function is approximated by: 

/^(x)^-^ (x»l/4). (3.74) 

V2 TT X 



That establishes an upper limit on Ad, 



A 



V2V^:>'-f, (3.75) 



where this approximation works. 



3.6 Parameters space. 

The model depends on 5 parameter. The medium height of PZ (L), the magnetic 
diffusion constant (Kq), the diffusion coefficient power index (5), the GW speed (K) 
and the Alfven velocity {Va) are determinated through the study of Nuclei CR measure- 
ments. Figure 13.81 and 13.91 show the results presented by Maurin et al. [Hj based on a 
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Maurin et al. (2001) 
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Figure 3.8: contour iso- lines for some of the parameters of the TZPM . The uncertainty 
bands for each parameters were obtained from the study of nuclei CR signals, spe- 
cially in the B/C ratio case ^4]. The two small plots show the contours when the 
power index 5 is fixed to 0.6. The big plot shows the contours when L = ?> kpc [64] . 

analysis to nuclei CR fluxes. To obtain that, they used the full propagation model, 
i.e. with all processes included, and solving this by means of sofisticated numerical methods. 

In the same spirit, Donato et al. [77] studied the problem of antiproton production 
obtaining constraints to the space of parameter (Table [3~T|) . which are totally compatible 
with B/C analysis. 

The case of PETE is simpler, in the sense, it depends on just three parameters - L, Kq 
and 5. Nevertheless, this reduced space of parameters is also described in the Maurin et al. 
framework because the leading processes, related to this three parameters, are taken into 
account for the positron and electron propagation. 



3.7 Solar Modulation 

When CR are closer to the Solar System, the Solar wind starts to pushing them away 
from the Solar neighborhood. The intensity in the cosmic-ray repulsion depends directly 
on Solar activity. The repulsion produces that CR lose energy and makes harder that CR 
can reach the Earth [75] . 
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A simple approach is the Force-Field Modulation [751 [HHl EI] , which models this phe- 
nomena as an electrical potential, which reduces CR energy and redistributes them. The 
CR intensity observed at Earth's orbit (Jg) is related with the interstellar intensity (Joo) 
as follows: 

2 2 

s — vn 

Jxie) = ^ ^ Joo(eoo), (3.76) 

where is the CR energy at the boundary of the heliosphere. 

Depending on the energy value, CR are affected differently: 

(3.77) 



V log ( — — — I + e + $ e < ffc 

p + e 



e + <& £c < £ 

where £c represent the limit where those are affected by two different modulation regime. 
The low-energy regime, e < Sc, is rarely used, however it is consistent with some studies 
about Solar Flares [71] . At a solar minimum, a typical value for is ~ 1 GeV [52] . 

The commonly used regime {Sc < s) depends on an energy shift which is related to the 
modulation potential (/>, 

$ = |Z|e0, (3.78) 

where \Z\e is the absolute value of CR electrical charge. The modulation potential, in 
principle, should vary on time. It can be estimated from diffusion coefficient and solar 



52 



3.7. SOLAR MODULATION 



Experiment 


year 


(f) [MV] 


MASS 89 


1989 


1400 ± 50 


MASS 91 


1991 


2000 ± 200 


CAPRICE 94 


1994 


664 ±5 


HEAT 94 


1994 


650 ± 50 


HEAT 95 


1995 


550 ± 50 


BETS 97+98 


1997, 1998 


600 ± 100 


AMS 98 


1998 


632 ± 13 



Table 3.2: Estimated values of solar modulation parameter for a selected set of CR experi- 
ments [85] . 

wind velocity. Although it is typically considered as a free parameter that should be fixed 
in each experiment |83j . 

Temporal correlations between Neutron Monitors (NM) measurements, the Solar ac- 
tivity and intensity of cosmic-rays have been found 84J. The information obtained in this 
type of analysis helps to improve the accuracy in the modulation potential (j) determination. 
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Chapter 4 



Positrons from DM annihilation 
in the galactic halo 

The production of positrons from DM annihilation is a very exciting pos- 
sibiUty to look for galactic DM. In order to study the positron signal it 
is necessary to study the propagation of positrons and the astrophysical 
uncertainties related to the propagation modeling. 

This chapter is based on our work [27] . 

4.1 Overview 

Secondary positrons and electrons are produced in the Galaxy from the interaction of 
nuclei CR on the ISM [39] and are an important tool for the comprehension of cosmic-ray 
propagation. Data on the cosmic positron flux (often reported in terms of the positron 
fraction) have been collected by several experiments [30l [85l [36l [371 EH [32] ■ 

The HEAT data pJIP have mildly indicated a possible excess of the positron fraction 
for energies above 10 GeV with respect to the available calculations for the secondary com- 
ponent [3H]. But recently, in October 2008, the lastest results of PAMELA experiment pS] 
have confirmed this feature [41j in the fraction. 

Different astrophysical contributions to the positron fraction in the 10 GeV region have 
been explored |30j , but only accurate and energy extended data could confirm the presence 
of a bump in the positron fraction and its physical interpretation. Alternatively, it has been 
conjectured that the possible excess of positrons found in the HEAT data could be due to 
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Halo model 


a 




7 


rs kpc 


Cored isothermal [SUj 


2 


2 





5 


Navarro, Frank & White [H] 


1 


3 


1 


20 


Moore [92] 


1.5 


3 


1.3 


30 



Table 4.1: Dark matter distribution profiles in the Milky Way. 



the presence of DM annihilation in the galactic halo [101 [HS]- Although, this interpretation 
is limited by the uncertainties in the halo structure and in the cosmic ray propagation 
modeling. 

In this chapter, the propagation of the positrons related DM annihilations are inspected 
through the solutions of the TZPM in connection with the study of the uncertainties due 
to propagation models compatibles with B/C measurements [Mi- To do this, the halo 
functions (Chapter H]) for some of the standard DM distributions have been calculated, as 
well, typical spectra of DM annihilation in beyond the Standard Model theories have been 
used to study the effects of the propagation uncertainties on the positron signal. 



4.2 Halo function and DM source 

We are here interested in primary positrons, namely the ones that are produced by the 
pair annihilations of DM particles. According to the various supersymmetric theories, the 
annihilation of a DM pair leads either to the direct creation of an electron-positron pair or 
to the production of many species subsequently decaying into photons, neutrinos, hadrons 
and positrons. We have considered four possible annihilation channels which appear in any 
model of weakly interacting massive particles (WIMP). 

The first one is the direct production of a e+e^ pair and is actually generic for theories 
with extra-dimensions like the UED models [571 IHHl [Ml ■ The energy of the positron line 
corresponds to the mass of the DM species. We have alternatively considered annihilations 
into W^W^ , T^T^ and bb pairs. These unstable particles decay and produce showers 
which may contain positrons with a continuous energy spectrum. 

Whichever the annihilation channel, the source term is generically written as 

.(x,s..) =,(..) (p(^y (4.1) 

where the coefficient 77 is a quantum term which depends on the particle being or not 
self-conjugate : for instance, for a fermion it equals 1/2 or 1/4 depending on whether 



56 



4.2. HALO FUNCTION AND DM SOURCE 



the WIMP is a self-conjugate particle or not. In what follows, we have considered a self- 
conjugate DM candidate (like the neutralino is supersymmetric theories) and taken t] = 
1/2. The annihilation cross section is averaged over the momenta of the incoming DM 
particles to yield (av) , the value of which depends on the specific SUSY model and is 
constrained by cosmology We have actually taken here a benchmark value of 2.1 x 
10^^^ cm'' sec^^ which leads to a relic abundance of fl^h^ ^ 0.14 (in agreement with the 
WMAP observations [94j[3]) under the hypothesis of dominant s-wave annihilation and by 
means of the relation: 

nh^ . 8.5 . 10-" = ^-^"'7-^-'^ (4.2) 

g*s[Xf) Xf-i{av) {av) 

where Xf = m^/Tf ~ (20-;- 25) with Tf the freeze-out temperature and where g-k{xf) 
and gi,s{xf) denote the effective number of degrees of freedom of the energy and entropy 
density at freeze-out, respectively. 

The DM mass is unknown. In the case of ncutralinos, theoretical arguments as well 
as the LEP and WMAP results constrain this mass to range from a few GeV [211 [Ml [S3 [H] 
up to a few TeV. Keeping in mind the positron HEAT excess, we have chosen a neutralino 
mass of 100 GeV. We have also analyzed the positron signal yielded by a significantly 
heavier DM particle of 500 GeV. 

Finally, all of the details in the energy distribution of the positrons ( -; ) produced 

in a single WIMP annihilation is described in the second part of Chapter [2l 

The only astronomical ingredient in the source term (Equation 14. l|) is the DM distribu- 
tion p(x) inside the Milky Way halo. We have considered the generic profile 

/rQX7 / l + (rQ/r,)" y^-^^/° 

where tq = 8.5 kpc is the galactocentric distance of the solar system. Notice that r 
denotes here the radius in spherical coordinates. The solar neighborhood DM density has 
been set equal to pQ = 0.3 GeV cm"'' [5^. We have discussed three profiles: an isothermal 
cored distribution [90] for which is the radius of the central core, the Navarro, Frenk and 
White profile [HI] (hereafter NEW) and Moore's model The NEW and Moore profiles 
have been numerically established thanks to N-body simulations. In the case of the Moore 
profile, the index 7 lies between 1 and 1.5 and we have chosen a value of 1.3 (Table |4~T]) . 



The possible presence of DM substructures inside those smooth distributions enhances 
the annihilation signals by the so-called boost factor, whose value is still open to de- 
bate [ini llOOj llOlj . Nevertheless, it has recently been shown that the boost factor due 
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Diffusion length A„ = (4K„t)i/2 [kpcj Diffusion length A„=(4KoT)1''2 [kpcj 

Figure 4.1: The halo convolution I is plotted as a function of the diffusion length Ad for various 
values of the PZ half-thickness L. The left panel features the case of an isothermal DM 
distribution whereas a NFW profile has been assumed in the right panel (Table [4Tl|) . 
When L is large enough for the positron horizon to reach the galactic center and its 
denser DM distribution, a maximum appears in the curves for Ad ~ rQ. 



to substructures in the DM halo depends on the positron energy and on the statistical 
properties of the DM distribution 102 . In addition, it has been pointed out that its 
numerical values is quite modest [103j . being of the order of 10-30. 

The positron flux at the Solar System is expressed as 

where the information related to particle physics is contained in 

«-77M (—) , (4.5) 

and the halo functions (/) for the different DM distribution are calculated using the 
methods described in Chapter (2] 



Figure ITTT] shows the effect of the overdensity at GC on the halo function. The two DM 
profiles are specially different in the central region, that becomes evident in the shape of 
each HF. When Ad is around 8.5 kpc - i.e. the distance to the galactic center - both HF 
present an increment that depends on how big is their respective DM profile in the GC 
surroundings. Furthermore, the value of L takes a very important role in the HF behaviors. 
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When L is small, both HF are similar because particles from the central region may escape 
earlier out of the PZ. On the other hand, a bigger L increase the chances of particles to 
arrive until the Solar System. 

However, we insist again on the fact that the true argument of the HF, whatever 
the approach followed to derive it, is the diffusion length Ad. The HF encodes the 
information relevant to cosmic ray propagation through the height L of the diffusive zone, 
the normalization A'o of the diffusion coefficient and its spectral index 5. It is also the only 
relevant quantity concerning the DM distribution. The analysis of the various astrophysical 
uncertainties that may affect the positron signal of annihilating WIMPs will therefore be 
achieved by studying the behavior of /. 



The central divergence 

Numerically derived DM profiles - NFW and Moore - exhibit a divergence at the center of 
the Milky Way. The density increases like r ' for small radii but cannot exceed the critical 
value for which the WIMP annihilation timescale is comparable to the age of the galactic 
bulge. The saturation of the density typically occurs within a sphere of ~ 10~^ pc, a much 
shorter distance than the size of the space increment in the numerical simulations. 

Fortunately, this numerical difficulty can be eluded by noticing that the positrons Green 
propagator does not vary much over the central DM distribution. This led us to replace 
inside a sphere of radius tq the cusp with the smoother profile 




where jo denotes spherical Bessel function of first kind 76J. The coefficients po, oi and 02 
are obtained by requiring that both the smooth density p* and its first derivative dp* /dr 
are continuous at ro. The other crucial condition is the conservation of the total number 
of annihilations within tq as the diverging cusp p on r ' is replaced by the distribution p* . 
These conditions imply that po = p{ro) whereas 

ai = 02 + 27, (4.7) 

" - 

In Figure 14.21 the halo integral / is plotted as a function of the diffusion length Ad in 
the case of the Moore profile and assuming a PZ half-thickness L ~ 10 kpc. Within the 
radius ro, the central DM divergence has been replaced either by a plateau with constant 
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Moore DM halo profile (r^ = 30 kpc) 




Diffusion length A„=(4Kot)i/2 [kpc] 



Figure 4.2: Same plot as Figure l4. 1[ the central DM profile within a radius ro is either a plateau 
at constant density po or the smooth distribution p* of Equation 14. 61 In the former 
case, the bump which / exhibits is significantly underestimated even for values of ro 
as small as 100 pc - solid dark blue - and drops as larger values are considered - 
solid light blue. On the contrary, if the DM cusp is replaced by the smooth profile p* , 
the halo integral no longer depends on the renormalization radius ro and the solid 
red and long-dashed black curves are superimposed on each other. 

density po or by the renormalized profile (Equation I4.6|l . In the case of the plateau, the 
maximum which / reaches for a diffusion length Ad ^ 7 kpc is underestimated even 
if values as small as 100 pc are assumed for tq. The larger that radius, the fewer the 
annihilations taking place within rg as compared to the Moore cusp and the worse the 
miscalculation of the halo integral. Getting the correct result featured by the solid red line 
would require a plateau radius so small that it would need to much numerical precision. 
On the contrary, we observe that the halo integral I is stable with respect to a change of 
ro when the renormalized density p* is used. 



4.3 Propagation uncertainties on the halo integral 

To compute the uncertainties in the positron signal, the best strategy is to compute 
first the HF for each DM profile, and in that way, we are equipped with a rapid enough 
method for scanning the ~ 1,600 different cosmic ray propagation models that have been 
found compatible [S3] with the B/C measurements. Each model is characterized by the 
half-thickness L of the diffusion zone and by the normalization Kq and spectral index d of 



60 



4.3. PROPAGATION UNCERTAINTIES ON THE HALO INTEGRAL 




Injection energy Eg [GeV] Injection energy Eg [GeV] Injection energy Eg [GeV] 



Figure 4.3: In each panel, the halo integral I is plotted as a function of the positron injection 
energy Ss whereas the energy e at the Earth is fixed at 10 GeV. The galactic DM 
halo profiles of Table 14. II are featured. The curves labeled as MED correspond to 
the choice of cosmic ray propagation parameters which best-fit the B/C ratio |64) . 
The MAX and MIN configurations correspond to the cases which were identified to 
produce the maximal and minimal DM antiproton fluxes [77], while the entire colored 
band corresponds to the complete set of propagation models compatible with the B/C 
analysis [64]. The MAX, MED and MIN configuration are given in Table [STTI 



the space diffusion coefficient. A large variation in these parameters is found in |64j and 
yet they all lead to the same B/C ratio. The height L of the PZ lies in the range from 1 
to 15 kpc. Values of the spectral index S extend from 0.46 to 0.85 whereas the ratio Kq/L 
varies from 10^^ to 8 x 10^^ kpc Myr^^ (See Chapter [3|). 

In the work about primary antiprotons performed by Donato et al. [77' , they found that 
the space of parameter obtained from the antiproton flux analysis is still compatible with 
the B/C analysis. Nevertheless, it is almost impossible to know if the same space produces 
equivalent results in the positron signal. In principle, nuclei CR and positrons (and 
electrons) are very different particles from the point of view of the physics of propagation. 
A first work on this topic was done by Hooper et al. [5^, but they scanned the TZPM 
parameters independently without any connection with other CR species analysis. 

In Figure 14.31 we have set the positron detection energy e at 10 GeV and varied the 
injection energy Eg from 10 GeV up to 1 TeV. The three panels correspond to the DM 
halo profiles of Table 14. II For each value of the injection energy e^, we have performed 
a complete scan over the 1,600 different configurations mentioned above and have found 
the maximal and minimal values of the HF / with the corresponding sets of propagation 
parameters. In each panel, the resulting uncertainty band corresponds to the yellow region 
extending between the two solid red lines. The lighter yellow domain is demarcated by 
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T. DelahQye, R. Lineros, N. Fornengo, F Donato & P. Saiatl (S007) 

Moore DM halo profile (r^, = 30 kpc) " 




10 100 1000 



Injection energy Eg [GeV] 

Figure 4.4: Same plot as Figure 14.31 the Moore DM profile has been selected. Four values of 
the positron detection energy e have been assumed. The flag-like structure of this 
figure results from the widening of the uncertainty band as the detection energy e is 
decreased. 

the long-dashed black curves labeled MIN and MAX and has a smaller spread. The MED 
configuration is featured by the long-dashed blue line. 

In Figure 14.41 the Moore profile has been chosen with four different values of the 
detection energy e. The corresponding uncertainty bands are coded with different colors 
and encompass each other as e increases. 

As Eg gets close to e, we observe that each uncertainty domain shrinks. In that regime, 
the diffusion length Ad is very small and the positron horizon probes only the solar 
neighborhood where the DM density is given by pq. Hence the flag-like structure of 
Figure 14.41 and a halo integral / of order unity whatever the propagation model. 

As is clear from Figure 14. H a small half-thickness L of the PZ combined with a large 
diffusion length Ad implies a small positron HF /. The lower boundaries of the various 
uncertainty bands in Figure l473l and |4 . 41 correspond therefore to parameter sets with i = 1 
kpc. Large values of Ad are obtained when both the normalization Kq and the spectral 
index 6 are large (Equation I3.28|) . However both conditions cannot be satisfied together 
once the B/C constraints are applied. For a large normalization Kq, only small values 
of 5 are allowed and vice versa. For small values of the detection energy e, the spectral 
index S has little influence on Ad and the conflguration which minimizes the halo integral 
/ corresponds to the large normalization Kq — 5.95 x 10~^ kpc^ Myr~^ and the rather 
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Lineros' PhD thesis (2008] 




Injection energy Eg [GeV] 



Figure 4.5: Diffusion length Ad versus injection energy. The Une at 5 kpc has been traced to 
denote that Ad bigger than this limit represent the range where the central DM 
overdensities start to appear. The two extreme B/C analysis configurations and the 
best-fit MED (Table [3 . II) have been used to calculated their associated Ad , showing 
how the MAX configuration covers longer distances instead of MED and MIN. 

small 6 ~ 0.55. For large values of e, the spectral index b becomes more important than 
if in the control of Ad ■ That is why in Figure 14.41 the lower bound of the red uncertainty 
domain corresponds now to the small normalization — 1.65 x 10^'^ kpc^ Myr^^ and 
the large spectral index b — 0.85. Notice that this set of parameters is very close to the 
MIN configuration of Table [3~T1 For intermediate values of e, the situation becomes more 
complex. Wc find in particular that for e = 30 GeV, the halo integral / is minimal for 
the former set of parameters as long as < 200 GeV and for the later set as soon as 
£s > 230 GeV. In between, a third propagation model comes into play with the intermediate 
values Kq = 2.55 x 10~^ kpc^ Myr~^ and 5 — 0.75. It is not possible therefore to single out 
one particular combination of Kq and 5 which would lead to the minimal value of the HF 
and of the positron DM signal. The MIN configuration which appeared in the antiproton 
analysis has no equivalent for positrons. 

The same conclusion holds, even more strongly, in the case of the upper boundaries of 
the uncertainty bands. Whatever the DM halo profile, the panels of Figure 14. II feature a 
peak in the halo function / for large values of L and for a specific diffusion length A^'^'^ ~ 7 
kpc. At fixed e and Eg, we anticipate that the maximal value for / will be reached for 
L = 15 kpc and for a diffusion length Ad as close as possible to the peak value X^'^^. Two 
regimes can be considered at this stage: 
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(i) When Ad < Xf^^^. This condition is satisfied when e and are close enough. 
Remember that e — Eg imphes Ad = for any set of Kq and S (Equation I3.28|) . 
Just to clarify the following discussion, Figure 14.51 shows the evolution of Ad for the 
parameter sets MAX, MED and MIN of the B/C analysis. The objective is to look 
for the largest possible values of Ad that maximize /. This is fulfilled at least by two 
propagation models. For small e, the large normalization Kq — 7.65 x 10~^ kpc^ 
Myr~^ is preferred with S = 0.46. We recognize the MAX configuration of Table [3~T] 
and understand why the long-dashed black curves labeled MAX in the panels of 
Figure 14.31 are superimposed on the solid red upper boundaries. For large e, the 
spectral index S dominates the diffusion length Ad and takes over the normalization 
Kq of the diffusion coefficient. The best model which maximizes / becomes then 
6 = 0.75 and Kq = 2.175 x 10"^ kpc^ nyT~\ 

(ii) When e and Eg are far enough, the diffusion length Ad may reach the critical value 
^max^ from values higher than Ag'^'^, for at least one propagation model which 
therefore maximizes the halo integral. As e and Eg are varied, the peak value of / is 
always reached when a scan through the space of parameters is performed. This peak 
value corresponds to the maximum of the halo integral, hence a horizontal upper 
boundary for each of the uncertainty bands of Figure [4~3] and [4~4l The set that leads 
to Ad = Xf^^^ is different for each combination of e and Eg and is not unique. In the 
case of the NFW DM profile of Figure 14.31 the halo integral / is maximized by more 
than 30 models above Eg = 120 GeV. 

The complexity of this analysis confirms that the propagation configurations selected 
by B/C do not play the same role for primary antiprotons and positrons. The two 
species experience the propagation phenomena, and in particular energy losses, with 
different intensities. As pointed out in Maurin et al. [104j . the average distance traveled by 
a positron is sensibly lower than the one experienced by an antiproton produced in the halo. 

4.4 Positron fluxes 

In the previous section we discussed the solutions of PETE and the astrophysical uncer- 
tainties on the HE. At this point, the next step is to apply this analysis on the positron 
signal, flux and fraction. As previously said, we do not adopt specific DM candidates, but 
do instead discuss the signals arising from a DM particle which annihilates into a pure 
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Model 


6 


Ko [kpcVMyr] 


L [kpc] 


MED 


0.70 


0.0112 


4 


Ml 


0.46 


0.0765 


15 


M2 


0.55 


0.00595 


1 



Table 4.2: Typical combinations of cosmic ray propagation parameters that are compatible with 
the B/C analysis [Bl]. The model MED has been borrowed from Table |3TT] Models 
Ml and M2 respectively maximize and minimize the positron flux for energies above 
10 GeV - the precise extent of which depends on the mass of the DM particle, on the 
annihilation channel and also on the DM profile. Note that Ml is the same as MAX 
in Table IsTT] but this is coincidental. 



final state. To do that, we consider four different specific DM annihilation channels: direct 
e^e^ production as well as W'^W^ , bb and r+r^. The DM annihilation cross section is 
fixed at the value 2.1 x 10"^^ cm'^ sec~^ and we consider the cases of a DM species with 
mass of 100 GeV and of 500 GeV. 

Generic DM candidates, for instance a neutralino or a sneutrino in supersymmetric 
models, or the lightest Kaluza-Klein particle in models with extra-dimensions, will entail 
annihilation processes with specific branching ratios into one or more of these benchmark 
cases. The positron flux in these more general situations would simply be a superposition 
of the results for each specific annihilation channel, weighted by the relevant branching 
ratios and normalized by the actual annihilation cross section. 

In Figure 14. 6| the propagated positron flux $e+ (times the positron energy squared) 
is featured as a function of e for a 100 GeV DM particle and a NFW density profile. 
The colored [yellow] area corresponds to the total uncertainty band arising from positron 
propagation. In all panels, the uncertainties increases at low positron energy. This may be 
understood as a consequence of the behavior of the HF / which was analyzed previously. 
Since positrons produced at energy and detected at energy s come from a sphere 
of radius Ad. That positron sphere grows as e decreases and so does the uncertainty 
band. On the other hand, the details of galactic propagation become less important in 
the fiux determination when positrons are originated closer to the Solar System. This is 
depicted in the HF because it goes to unity independently from the propagation parameters. 

Notice also that the uncertainty band can be sizeable and depends significantly on the 

positron spectrum at production — . In the case of the e+e^ line of the upper left 

panel, the positron flux (i>g+ exhibits a strongly increasing uncertainty as e is decreased 
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Figure 4.6: Positron flux E^$^+ versus the positron energy e, for a DM particle with a mass of 
100 GeV and for a NFW profile (Table |4TT|) . The four panels refer to different annihi- 
lation final states : direct e'^e~ production (top left), bb (top right), W^W~^ (bottom 
left) and t'^t~ (bottom right). In each panel, the thick solid [red] curve refers to the 
best-fit choice (MED) of the astrophysical parameters. The upper [blue] and lower 
[green] thin solid lines correspond respectively to the astrophysical configurations 
which provide here the maximal (Ml) and minimal (M2) flux - though only for en- 
ergies above a few GeV in the case of (Ml). The colored [yellow] area features the 
total uncertainty band arising from positron propagation. 
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Figure 4.7: Same plot as in Figure 14 . 61 but with a DM particle mass of 500 GeV. 
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from niy^ down to 1 GeV. That uncertainty is one order of magnitude at e = 10 GeV, and 
becomes larger than 2 orders of magnitude below 1 GeV. Once again, the positron sphere 
argument may help to understand the situation. At fixed detected energy e, the radius Ad 
increases with the injected energy Sg- We therefore anticipate a wider uncertainty band as 
the source spectrum gets harder. This trend is clearly present in the panels of Figure 14.61 
Actually direct production is affected by the largest uncertainty, followed by the r+r^ 
and W^W~ channels where a positron is produced either directly from the or from 
the leptonic decays. In the hh case, which is here representative of all quark channels, a 
softer spectrum is produced since positrons arise mostly from the decays of charged pions 
originating from the quark hadronization. Most of the positrons have already a low energy 
Eg at injection and since they are detected at an energy e ~ e^, they tend to have been 
produced not too far from the Earth, obtaining smaller dependence to the propagation 
uncertainties. 

The astrophysical configuration M2 (Table 14. 2p provides the minimal positron fiux. It 
corresponds to the lower boundaries of the yellow uncertainty bands of Figure 14.61 The 
Ml configuration maximizes the flux at high energies. For direct production and to a lesser 
extent for the r+r~ channel, that configuration does not reproduce the upper envelope 
of the uncertainty band in the low energy tail of the flux. As discussed in the previous 
section, the response of $£+ to the propagation parameters depends on the detected energy 
£ in such a way that the maximal value cannot be reached for a single astrophysical 
configuration. Finally, taking as a reference the median fiux, the uncertainty bands extend 
more towards small values of the fiux. In all channels, the maximal fiux is typically a factor 
of ~ 1.5-2 times larger than the median prediction. The minimal fiux features larger devia- 
tions with a factor of 5 for the hh channel at E' = 1 GeV, of 10 for W^W~ and of 30 for r+r" . 

Figure 14.71 is similar to Figure 14.61 but with a heavier DM species of 500 GeV 
instead of 100 GeV. Since the mass is larger, on average the injected energy is 
higher, too. Notice that at fixed positron energy e, the radius Ad of the positron sphere 
increases with Sg- And then, we anticipate that the propagated fluxes are affected by 
larger uncertainties for heavy DM particles. Again, the maximal fiux does not exceed 
twice the median flux, while the minimal configurations are significantly depressed. 
At the reference energy e = 1 GeV, reductions by a factor of 10 between the median 
and minimal predictions are obtained for the hh channel and amount to a factor of 20 in 
the W^W^ case. They reach up to 2 orders of magnitude for the direct positron production. 

In this large DM mass regime, the astrophysical configuration M2 does not reproduce 
by far the lower bound of the uncertainty band as it did for the 100 GeV case. That draws 
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two possibilities. 

(i) Once the full positron spectrum (not each single channel) at the source is chosen, the 
correct determination of the uncertainty in the flux requires a full scan of the prop- 
agation parameter space for each energy e. The use of representative astrophysical 
configurations such as Ml and M2 would not provide the correct uncertainty over 
the entire range of positron energy e . 

(ii) Specific predictions have to be performed for a given model of DM particle and 
a fixed set of astrophysical parameters. This is why fits to the experimental data 
should be performed for each propagation configuration over the entire range of 
the measured positron energies e. The best fit should correspond to a unique set of 
astrophysical parameters. This procedure is the only way to reproduce properly the 
correct and specific spectral shape of the flux. 

The effect induced by different DM proflles is presented in Figure 14.81 where the 
positron fluxes for the bb and W^W~ channels are reproduced for the DM distributions 
of Table 14. 11 The mass of the DM particle is fixed at = 100 GeV. Notice how steeper 
profiles entail larger uncertainties, especially for the upper bound. This is mostly due to 
the fact that for large values of L the positron flux is more sensitive to the central region 
of the Galaxy, where singular profiles like the NFW and Moore distributions have larger 
densities and therefore induce larger annihilation rates. On the contrary, the lower envelope 
of the uncertainty band is not affected by the variation of the halo profile. In this case, 
with typically small heights L, positrons reach the solar system from closer regions, where 
the three halo distributions are very similar and do not allow to probe the central part of 
the Milky Way. 

Figure 14.91 depicts the information on the positron fiux uncertainty from a different 
perspective. The fiux <i>g+ and its uncertainty band are now featured for fixed values of the 
detected energy e whereas the DM particle mass is now varied. The flux $g+ is actually 
rescaled by the product <i>g+ for visual convenience. Each band corresponds to a 

specific detected energy e and consequently starts at = e. 

In the case of the W~^W~ channel, the bands start at = mw because this channel 
is closed for DM masses below that threshold. The behavior of these bands can be 
understood from Figure 14.41 where the halo function / is plotted for the same detected 
energies, as a function of the injection energy eg. In the case of direct positron production, 
there is a simple link between the two figures, because the source spectrum in this case is 
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Figure 4.8: Positron flux e^^^+ versus the positron energy e, for a DM particle mass of 100 GeV 
and for different halo density profiles : cored isothermal sphere [90] (left panels), 
NFW i91j (central panels) and Moore [52] (right panels) (Table 147111 ■ The upper 
and lower rows correspond respectively to a bb and W^W^ annihilation channel. 
In each panel, the thick solid [red] curve refers to the best-fit choice (MED) of the 
astrophysical parameters. The upper [blue] and lower [green] thin solid lines stand 
for the astrophysical configurations Ml and M2 of Table 14 . 21 The colored [yellow] 
area indicates the total uncertainty band arising from positron propagation. 
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Figure 4.9: For fixed values of the detected energy e, the uncertainty bands on the positron flux 
£^m^$j,+ are shown as a function of the mass of the DM particle. The energies 
considered in the flgure are £ = 3, 10, 30 and 100 GeV. Each band refers to one of 
those values and starts at = e . 
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just a line at Eg = rn^. For the other channels the situation is more involved since we 
have a continuous injection spectrum with specific features as discussed above. The main 
information which can be obtained from Figure 14.91 is that at fixed detection energy, the 
larger the DM mass, the larger the uncertainty. Let us take for instance a detection energy 
of ii' = 3 GeV. For direct production, where Ss = m^, increasing the DM mass translates 
into a larger radius Ad of the positron sphere. As a consequence, the uncertainty band 
enlarges for increasing masses. This occurs for all the annihilation channels, but is less 
pronounced for soft spectra as in the bb case. Similar conclusions hold for all the other 
values of e. 



In addition, to compute a full positron signal, we need to deal with the secondary 
positrons and the full electron signal. To simplify the situation we used the results of 
Moskalenko et at. [39], which were parameterized by Baltz et al. [40] obtaining function for 
these fluxes: 



d$\ 0.16 e- 



de 1 + 11 + 3.2 £2.15 

d^y 0.7 £0-7 



de ) 1 + 110 £1-5 + 600 £2-9 + 580 £" 

d^Y^ 4.5 £0-7 



de 1 + 650 £2-3 + 1500 £4-2 



(4.9) 
(4.10) 
(4.11) 



where £ = e/1 GeV and all those are expressed in units of GeV ^ cm 2 sec ^ sr ^. Notice 
that the full electron signal is composed by the addition of secondaries and primaries, 

U'^-X . (4.12) 



Comparison with available data is presented in Figure 14. 101 14. Ill and 14. 121 In Fig- 
ure 14. 10| the positron fraction 



$TpT 



$TOT , ^TOT 
e~ e+ 



(4.13) 



is plotted as a function of the positron energy e . The total positron flux $ J-P""" at the Earth 
encompasses the annihilation signal and a background component ~ see the thin solid 
[brown] lines. The mass of the DM particle is 100 GeV and a NFW profile has been assumed. 
The data from HEAT [3D], AMS [31137], CAPRICE [31] and MASS [31 are indications of 
a possible excess of the positron fraction for energies above 10 GeV. Those measurements 



72 



4.4. POSITRON FLUXES 



— I — I — I I 1 1 1 1 1 1 — I — I I 1 1 1 1 1 

NFW Halo profile (r, = 20 kpc) 

<av> = 2.1 X 10""' cm"^ s~' 
: 100 GeV 



T. Delahaye. R Lineros. N. Form 
1 1 ^1 II I I ■ f 



I, F. Donato & P.Salati |2007) 



-I 1 I I I I I I 





Direct prod. 

Boost factor = 10 




B/C best fit 
uncer. band 
background 



bb channel 

Boost factor = 50 



I I I ll| 




Heat 2000 
MASS-91 
CAPRICE94 



W"'"W channel 

Boost factor = 30 



_l 



-I 




z\ channel 

Boost factor = 40 
[I I I I 



10" 



10' 

Positron energy [GeV] 



10" 



10" 



10' 

Positron energy [GeV] 



10" 



Figure 4. 10 : Positron fraction e+/(e^ +e"^) versus the positron detection energy e. Notations are 
as in Figure [476] In each panel, the thin [brown] solid line stands for the background 
|40l |39] whereas the thick solid [red] curve refers to the total positron flux where the 
signal is calculated with the best-fit choice (MED) of the astrophysical parameters. 
Experimental data from HEAT l20j, AMS [31 [37], CAPRICE [ST] and MASS [32] 
are also plotted. 
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Figure 4.11: Same plot as in Figure [4 . 101 but with a mass of the DM particle of 500 GeV. 
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may be compared to the thick solid [red] line that corresponds to the MED configuration. 
In order to get a reasonable agreement between our results and the observations, the 
annihilation signal has been boosted by an energy-independent factor ranging from 10 to 
50 as indicated in each panel. At the same time, the positron background - for which we 
do not have an error estimate yet - has been shifted upwards from its reference value of 
Baltz et al. (Equation 14 .911 by a small amount of 10%. 

As is clear in the upper left panel, the case of direct production offers a very good 
agreement with the potential HEAT excess. Notice how well all the data points lie within 
the uncertainty band. A boost factor of 10 is enough to obtain an excellent agreement 
between the measurements and the median flux. A smaller value would be required for 
a flux at the upper envelope of the uncertainty band. The W^W^ and r+r^ channels 
may also reproduce reasonably well the observations, especially once the uncertainty is 
taken into account, but they need larger boost factors of the order of 30 to 40. On the 
contrary, softer production channels, like the bb case, are unable to match the features 
of the putative HEAT excess for this value of the DM particle mass. For all annihilation 
channels, the uncertainty bands get thinner at high energies for reasons explained above. 
They surprisingly tend to shrink also at low energies, a regime where the positron horizon 
is the furthest and where the details of galactic propagation are expected to be the most 
important. Actually, the annihilation signal turns out to be completely swamped in the 
positron background. In particular, the signal from direct production stands up over the 
background only for energies larger than 5 GeV. The corresponding uncertainty on the 
positron fraction is at most of the order of 50% for energies between 10 and 20 GeV. In the 
other cases, the uncertainty bands are even thinner. 

Note that the positron background also should be affected by uncertainties due to 
secondary production processes and propagation. This chapter is just devoted to study 
the uncertainties at the level of primary positron production. However, the study of 
uncertainties of the secondary production is in Chapter [s] 

Somehow different is the situation for larger masses of the DM candidate. Figure 14. Ill 
features the same information as Figure 14. 101 but now for = 500 GeV. In this case, 
all the annihilation channels manage to reproduce the experimental data, even the softest 
one bb. For direct production, the positron fraction is very large at energies above 40 GeV, 
where no data are currently available. This feature would be a very clear signature of 
DM annihilating directly into e+e" pairs, with strong implications also on the nature of 
the DM candidate. For instance, bosonic dark matter would be strongly preferred, since 
Majorana fermionic DM, like the neutralino, possesses a very depressed cross section into 
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Figure 4.12: Positron flux E^^^+ (not fraction) versus the positron energy e, for a 500 GeV 
DM particle. Notations are the same as in Figure 14 . 101 Experimental data from 
HEAT 130], AMS [MIEZ], CAPRICE [31] and MASS 132. are plotted. 
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light fermions because of helicity suppression in the non-relativistic regime. Astrophysical 
uncertainties on the signal in this case show up more clearly than for the case of a lighter 
DM species, but still they are not very large. The drawback of having a heavier relic is 
that now the boost factors required to match the data are quite large. In Figure 14. Ill 
they range from 250 for the soft channel to 400 for the r+r^ case. Such large boost 
factors appear to be disfavored, on the basis of the recent analysis of Lavalle et al. [10211103] . 

In Figure 14. 12| the positron flux (not the fraction) is compared to the available 
experimental data for a 500 GeV DM particle and a NFW profile. The solid thin [brown] 
line features the positron background which we shifted upwards by 10% with respect to 
the reference value of Baltz et al. [3D]. The thick solid [red] line encompasses both that 
background and the annihilation signal which we calculated with the best-fit choice (MED) 
of the astrophysical parameters. Both curves have been derived assuming solar modulation 
implemented through the force field approximation with a Fisk potential (pp of 500 MV 
(chapter [3|). The dashed [red] line instead corresponds to the total positron interstellar 
flux without solar modulation. Notice that this curve is superimposed on the thick [red] 
line above ~ 10 GeV, a regime where cosmic ray propagation is no longer affected by the 
solar wind. A reasonably good agreement between the theoretical predictions and the data 
is obtained, especially once the theoretical uncertainties on the annihilation signal are 
taken into account. Notice that the spread of each uncertainty band is fairly limited as we 
already pointed out for the positron fraction. The reasons are the same. 

Prospects for the future missions are shown in Figure 14. 131 and 14. 141 In Figure 14. 13| a 
100 GeV DM particle and a NFW halo profile have been assumed. The median [red] curve 
corresponds to the prediction for the best-fit MED choice of astrophysical parameters 
whereas the upper [blue] and lower [green] lines correspond respectively to the Ml and M2 
propagation models (Table 14 . 2|) . Since we are dealing with predictions which will eventually 
be compared to the measurements performed over an entire range of positron energies, 
we have to choose specific sets of propagation parameters. The upper and lower curves 
therefore do not represent the maximal uncertainty at each energy - though they may do 
so in some limited energy range - but instead they are "true" predictions for a specific 
set of propagation parameters. Figure 14. 131 summarizes our estimate of the capabilities 
of the PAMELA detector [38 after 3 years of running. We only plotted statistical errors. 
We reached the remarkable conclusion that not only will PAMELA have the capability 
to disentangle the signal from the background, but also to distinguish among different 
astrophysical models, especially for hard spectra. Our conclusion still holds for the bb soft 
spectrum for which the Ml, MED and M2 curves of the upper right panel differ one from 
each other by more than a few standard deviations. PAMELA could be able to select 
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among them, even when systematical errors are included. 

In Figure 14. 141 the case of a 500 GeV DM particle is confronted with the sensitiv- 
ity of AMS-02 for a 3-year flight. The possibility to disentangle the signal from the 
background is also clearly manifest here, even once the astrophysical uncertainties are 
included - provided though that boost factors of the order of 200 to 400 are possible. 
But, unless direct production is the dominant channel, a clear distinction among the 
various astrophysical models will be very difficult because the Ml and M2 configurations 
are closer to the MED curve now than in the previous case of a lighter DM species. 
Comparison between Figure 14. 131 and Figure 14. 141 clearly exhibits that at least below the 
TeV scale, the effect of the mass rriy^ should not limit the capability of disentangling the 
annihilation signal from the background. More problematic is our potential to distinguish 
among different astrophysical models when the DM mass sizeably exceeds the 100 GeV scale. 

Note added: As the time of our publication [57], on which this chapter is based, the 
PAMELA data of the positron fraction were not yet released. Nevertheless, most of the 
above analysis is still valid with the current state of the positron excess feature, which 
will be discussed in the next chapter in conection with the calculation of the positron 
background. 
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Figure 4.13: Predictions for PAMELA for a 3-year mission. The positron fraction e'^/(e~ + e"^) 
and its statistical uncertainty are plotted against the positron energy e for a 100 GeV 
DM particle and a NFW profile. Notations are the same as in Figure 14. 101 The 
thick solid curves refer respectively to the total positron flux calculated with the 
Ml (upper [blue]), MED (median [red]) and M2 (lower [green]) sets of propagation 
parameters. 
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Figure 4.14: Predictions for AMS-02 for a 3-year mission. The positron flux e^"l>e+ f^nd its 
statistical uncertainty are featured as a function of the positron energy e for a 
500 GeV DM species and a NFW profile. Notations are the same as in Figure 14. 121 
The thick solid curves refer respectively to the total positron flux calculated with the 
Ml (upper [blue]), MED (median [red]) and M2 (lower [green]) sets of propagation 
parameters. 
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Chapter 5 



Secondary positron flux at 
Earth 

Positrons in the galactic environment can be produced in many ways. 
The main component in the positron signal comes from nuclei cosmic-rays 
interactions with the interstellar medium. 

This chapter is partially based on our work [105] . 

5.1 Overview 

Cosmic positrons are created by interaction of cosmic-ray nuclei on interstellar matter and 
propagate in a diffusive mode, because of their interaction with the turbulent component 
of the galactic magnetic field. The expected flux of positrons can be calculated from 
the observed cosmic ray-nuclei fluxes, using the relevant nuclear physics and solving the 
transport equation. 

The HEAT experiment [30] showed that the positron fraction possibly exhibits an 
unexpected bump in the 10 GeV region of the spectrum. Although it is not excluded that 
this bump could be due to some unknown systematic effect, the HEAT result has triggered 
a lot of explanations. For instance, Moskalenko et al. [31] suggested that an interstellar 
nucleon spectrum harder than the expected one could explain the excess. A lot of works also 
focused on the dark matter hypothesis, the bump being due to a primary contribution from 
the annihilation of dark matter particles. The positron excess expected in this framework 
is very uncertain, because the nature of dark matter is not known, and because the prop- 
agation of positrons involves physical quantities that are also not currently precisely known. 
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The related astrophysical uncertainties were calculated and quantified in [27j . where 
it has been shown that they may be sizeable, especially in the low energy part of the 
spectrum, a property which is common also to the antiproton [77] and the antideuteron 
[ins [MI signals. 

For positrons, sizeable fluxes from dark matter annihilation are typically possible 
if dark matter overdensities are locally present, a fact that is usually coded into the 
so-called "boost factor". A detailed analysis on the admissible boost factors for positrons 
and antiprotons has been performed by Lavalle et al. [108j . who have shown that boost 
factors are typically confined to be less than about a factor of 10-20. Computing the 
antimatter fluxes directly in the frame of a cosmological N-body simulation leads to the 
same conclusions [109] . 

The PAMELA experiment |3H] has just released its first results on the positron fraction 
for energies ranging from 1.5 GeV to 100 GeV and with a large statistics jH]. The 
positron fraction is observed to steadily rise for energies above 10 GeV, reinforcing the 
possibility that an excess is actually present. The calculation of the uncertainties affecting 
the standard spallation-induced positron population is in fact especially important when 
unexpected distortions are observed in experimental data, in order to properly address 
the issue. This, together with the recent analysis on the positron signal from dark matter 
annihilation and its astrophysical uncertainties [27] . will set the proper basis to discuss in 
details the experimental results. 

The uncertainties on the positron flux have several origins: 

i) The cosmic-ray nuclei measurements come with their experimental uncertainties, 
which then affect the predictions of induced secondary fluxes, like positrons. 

ii) Various modelings of the nuclear cross sections involved in the positron production 
mechanism are available, and they do not exactly match each other, implying a range 
of theoretical variation. 

iii) The uncertainties in the propagation parameters involved in the transport equation 
have been thoroughly studied jM]. 

In this chapter, we study the production of secondary positrons which come from the 
interaction of proton and alpha particles cosmic-rays with the interstellar medium formed 
by Hydrogen and Helium. As well, the propagation of positrons is modeled according to the 
TZPM. Finally, we calculate and analyze the positron flux and fraction comparing them 
with available experimental data. In each part the uncertainties related to propagation. 
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nuclear cross section and cosmic-rays measurement were studied. 



5.2 Production of secondary positrons. 

Secondary positrons are produced by interactions of Nuclei CR, principally protons and 
alpha particles, with the ISM composed mainly by Hydrogen and Helium. 

The source term of PETE (Equation 13. 17|l is generically described by: 

Scab(x,ec) = 47r na$b;7^(a + b c + X) , (5.1) 

where ria is the number density of species a, <i>b is the total flux of species b per solid angle 

and — — represents the differential cross section for the specific process. 
de 

In the case of secondary positrons, the source term related to the scattering between 
protons and Hydrogen is more complex than the former one, because we are in presence of 
a proton flux which continuously bombards the ISM. In this case, the source term is: 

/A I \ 
dSp <I>p(x, Ep) (gp,ge+) , (5.2) 

where nn is the Hydrogen number density, {r]a) represents the inclusive cross section 
of the process p + p ^ e+ + X discussed and calculated in Chapter [2l Let us specify 
that the inclusive cross section, in our case, comes from three invariant cross section 
parameterizations Badhwar et al. [48j, Tan and Ng [5^ and Kamae et al. [55] . 

The proton flux (per solid angle) $p - in principle - depends on the position, although 
the effects of magnetic diffusion help to homogenize the CR density, in addition to that, 
positrons (electrons) just can propagate not so far from the production point due to the huge 
energy losses. For this reason, we can assume that the proton (alpha particles) CR density 
is homogeneous in all the PZ [105j . To avoid the problem of modeling the propagation of 
Nuclei CR in the Galaxy, we used flux parameterizations of proton and alpha particles 
based on experimental data. The Shikaze et al. parameterization [17] , 

$(e) = A (^^) sec-i sr-^ (GeV/n)-i) , (5.3) 

is a function of the rigidity and is based on BESS experiment results. On the other hand, 
the Donato et al. parameterization [65j . 

$(£) =N ( £Mii!£il^ , (5.4) 



1 GeV/n 
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Donato et al. [65] parameterization 




proton 


alpha 


TV 


1.3249 X lO'^ 


7.21 X 10^ 


7 


2.72 


2.74 


Shikaze et al. [47^ parameterization 




proton 


alpha 


A 


(1.94 ±0.13) X 10"* 


(7.10 ±0.56) X 10^ 


Pi 


0.70 ±0.52 


0.50 ±0.31 


P2 


2.76 ±0.03 


2.78 ±0.03 



Table 5.1: Values of parameters for Donato et al. and Shikaze et al. parameterizations of the 
proton and alpha particles fluxes. The units of fluxes for each parameterization are 
m"^ sec"^ sr"^ (GeV/n)"\ 



is a function of the kinetic energy per nucleon and is based on data from BESS and AMS 
experiments. Notice that the functional forms are quite different, but both describe well 
enough the available data. The values of parameter are given in Table 15 . II 

Figure 15. II shows the effects produced by different proton flux parameterizations. 
The curves calculated with the Shikaze's parameterization are qualitative equivalent 
to the ones calculated with the Donato's parameterization. On the other hand, most 
of the differences come from the nuclear cross section parameterizations. We see that 
Kamae's parameterization produces lower values of the source term instead of Tan 
and Ng and Badhwar ones, which have a lower relative error among them. At this 
level, all these differences can be considered like uncertainties which come from the nu- 
clear physics context, because all of those were produced from equivalent experimental data. 

The complete source term, which considers the interaction of protons, alpha particles 
with Hydrogen and Helium, has the same shape that the proton-Hydrogen case: 

Qfuu(x,ee+) = 47r ^ ^ ^^(x) /" de^ $j(£j) ^^y^(£j, £e+) , (5.5) 

where the inclusive cross sections for processes p+He, a+H and a+He have been described 
according to scaling factors s/, 

where Aj is the mass number of the incident particle. The scaling factors could be quite 
controversial in the nuclear physics field, however, those reproduce good enough the avail- 
able experimental data. Basically, we used two of the most known scaling factors. The Orth 
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Figure 5.1: Positron source term s^'^^q versus positron energy, 
et al. scaling factor [llQj . 

and the one proposed by Norbury et al. 



(5.7) 



(5.8) 



Numerically speaking, both scaling factors produce slightly different values. For example 
the Orth's scaling factor is 2.54 instead of 2.76 from Norbury's factor for the case of p+He. 
That difference is screened in the final result due to the low intensity of alpha particles - 
in comparison to proton flux - and the low density of Helium in the ISM. 



Another point that remains unseen is the spatial distribution of Hydrogen and Helium 
in the ISM. The interstellar gas principally is concentrated in the Galactic Plane. In first 
approximation, the gas is homogeneously distributed with densities njj = 0.9 cm^'^ and 
"-He = 0.1 cm^"^ jJB]. Many structures are also present in the GP like the spiral arms, the 
galactic bulge, among others. However, positrons and electrons come from local distances, 
where the gas density is more or less constant, for these reason the gas in the galactic plane 
is considered homogeneous in the thin disk of the TZPM (Chapter [3]), i.e. in a cylinder of 
height 2hz. In other words, the source term for secondary positron is the Equation 15 . 51 but 
with a gas density distribution given by: 



"-H/Hc(x) 



"H/Hc \z\ < 

other cases 



(5.9) 
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Figure 5.2: Halo function for the secondary positron source term versus the diffusion length Ad 
(left). Corrective function fc versus Xo/L^. Both the numerically calculated function 
and the analytical fit are plotted (right). 

where n-n/Hc correspond to the values previously given. 

5.3 Propagation of secondary positrons. 

As it was already discussed m the previous chapters, specially in Chapter |3l the propagation 
of positrons is modeled according to the TZPM. The HF associated to the secondary 
positron source term (Equation 15 . 5|l is based on the HF for homogeneous sources and is 
easier to calculate (Equation l3.68|) . In our case, R and i?si.c are set to 20 kpc and to 100 pc. 

Figure 15.21 show the halo function for the secondary positron source term. This HF 
presents a fast decreasing behavior for values of Ad bigger than /i^, which is related to the 
propagation volume and the amount of sources inside it. In other words, when Ad < 
positrons come from a propagation sphere with radius Ad which is full of sources. When 
Ad goes beyond this limit, the sphere starts to be partially full. 

Unlike the dark matter HF case, the effects of vertical boundaries play a less important 
role. Those still modify the HF for big values of Ad, enhancing its decreasing behavior. 
Nevertheless, this effect appears in the range where HF is small, expecting a smaller 
contribution in the flux (Figure 15 . 2|) . The three curves show how vertical boundaries speeds 
up the decreasing behaviors when Ad is bigger than L^- As well, all the curves have a 
common behavior until boundary effects appear. This feature, that they have in common. 
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Parameters fc 




ai a2 as as 


6o bi 62 


0.6 


0.0212 0.0267 -1 0.5 -0.0868 


0.9713 -0.5046 -0.6229 



Table 5.2: Parameter values for the fit of the fc function. 



can be used to produce an unique representation of them. 

The corrective function fc is a transfer function from the boundary-free HF respect to 
the one with vertical boundaries (right-panel in Figure 15 . 2|) , 

7scc.(AD,i2) = fci>^D/L^) X 7scc.(AD,i2 = oo) . (5.10) 

It appears to be common to any HF, in the specific case of secondary positrons. Also, this 
nice property helps to improve the speed in flux calculations, because the only function that 
should be computed is the boundary-free HF. The corrective function can be parameterized 
as follows (Figure I5.2|l , 

{1 X < Xr 

eyip {ait + a2t^ + a'it^ + a^t^ ^ a^t^) Xr<x<i , (5.11) 
exp (bo + bit + b2t^) i<x 

where t = x — Xr and the values of the parameters are given in Table 15 .21 

Up to this point, the secondary positron propagation has been discussed through 
the analisis of the halo function. The following step is to calculate the flux at the Solar 
System position. As it was seen in Equation 13. 581 the flux is obtained by the convolution in 
energy between the source term and the halo function. However this procedure should be 
performed with proper care due to the fast changes that the halo function has. According 
to Figure 15 . 1\ the source term can be roughly assumed as a power-law with spectral index 
-2.7. On the other hand, Figure [5721 shows that the HF is equal to 1 when Ad ^ 0.1 kpc 
and starts to decrease as Ad grows. This feature, when translated into energies, produces 
an extremely fast change in the function to be integrated, as Figure |4 . 51 shows. 

From the value of Ad and the observed energy e, it is possible to infer the value of the 
injected positron energy Cg-. 

, , 2 \ 1/(5-1) 

e.= (^e^-i + (<5-l)^j , (5.12) 

where A§ — KqTe- For example, using the MED set (Table 14.21) . a diffusion length of 
0.1 kpc and observed energies of 1 and 10 GeV, we obtain injection energies of 1.003 and 
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Figure 5.3: Interstellar electron and positron fluxes e <l>e for the MED set (Table 13 . l|l versus 
energy. The curves correspond to fluxes calculated from Nuclei CR interactions with 
the ISM. Each curve represents a different nuclear cross section and Nuclei CR flux 
parameterization. Also the Strong et al. (Equation 14 .Qp flux parameterizations are 
shown. The uncertainty band related to those parameterization is plotted (yellow 
band) as well. 



10.056 GeV respectively. That means the main contribution of the source term come from 
the range between the observed energy and the injected energy at 0.1 kpc. After that hmit, 
the halo function and the source term decreases fast, contributing less to the final results. 



5.4 Positron flux and uncertainties. 

Current experiments, as HEAT [30], AMS [37j among others, are able to measure the 
positron and electron fluxes that arrive to the Earth. As well, we are able to compute the 
positron and electron flux using the propagation model and sources that we discussed in 
last sections. 
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Figure [5~3] shows the positron and electron fluxes calculated using the two nuclei CR pa- 
rameterization (Donato et al. [HS] and Shikaze et al. [Uj) and the three nuclear cross section 
parameterizations (Tan and Ng [50], Kamae et al. [55] and Badhwar et al. [48]). The yel- 
low band represent the uncertainties from nuclear cross sections and nuclei CR observations. 

For the different calculated fluxes, we observed that at high energy (> 100 GeV) all fluxes 
converges into one common curve. This is produced by a common asymptotic behavior of 
all parameterizations. On the other hand, differences arise at low energies because smaller 
differences from high energies start to accumulate when the flux is calculated. Let us to 
remark that the secondary positron and electron fluxes depend on an integration on energy 
from inflnity down to the observed energy. 

Moreover, in Figure 15.31 the reference curves of Baltz et al. parameterization 
(Equation I4.9|l . which come from the results of Strong and Moskalenko [3H], are plotted. 
We note how positron calculated with our model exhibits same power-law behavior with 
power index around -3.5 like Strong et al. For energies below 10 GeV, bigger differences 
appear. However, that is the range of energies dominated by the Solar modulation, which 
indirectly reduces the differences among the modulated fluxes. The secondary electron 
flux presents bigger difference with respect to secondaries from Strong et al. However, the 
secondaries contribution to electron flux is highly diluted into the sea of primary electron 
flux. This establishes a first motivation for further improving calculations in other to 
estimate the primary electron flux, as well. 

Extra information is obtained from the study and analysis of the uncertainty in 
the propagation model. In a similar way to what was done in the case of dark matter 
annihilation, we calculate the positron flux considering the 1600 parameter configurations 
compatibles with the B/C measurement [64] . 

Figure 15.41 and Figure 15.51 show the positron fiux e'^-^<i>g+ calculated for the MED 
set and in the solar-modulated and -unmodulated regime (solid and dashed red lines, 
respectively), the reference curves for Strong et al. secondary positrons are plotted for the 
same regimes. 

In both figures, uncertainty bands associated to different propagation model are 
included. The upper and lower bounds in each band correspond to the maximum and the 
minimum fluxes obtained among all the parameter space. In a similar way that in the case 
of DM annihilation, there are not unique parameter sets which maximize or minimize the 
flux in the whole range of energies. 
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Figure 5.5: Secondary positron flux £^'^<E>g+ versus positron energy. Positron fluxes were calcu- 
lated using proton and alpha CR fluxes from Donato et al. [65] and Shikaze et al. [47] 
with the Tan et al. cross section parameterization [50J. 
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The first interesting issue that we notice is how the experimental data is nicely contained 
inside the uncertainty band, this fact corroborates the consistency of the propagation 
model among different species of cosmic rays. Furthermore, we note that the fluxes at 
low-energy (< 10 GeV) present similar values due to the solar modulation effects. 

From the figures, we see how different nuclei CR flux parameterizations do not modify 
too much the positron flux. For example, in Figure 15 .41 we compare the positron fluxes 
calculated with the same cross section parameterization, this case Kamae's one, but with 
each of the CR flux parameterizations, Donato's and Shikaze's one. The same situation 
occurs in Figure 15.51 

The major differences arise in the low-energy range but just for interstellar fluxes. These 
differences cannot be easily seen in the modulated signal because the Solar modulation 
generate an energy shift equivalent to 0.6 GeV respect to the energy of the interstellar 
flux. This energy shift is enough to allow us to see at Earth just the range of energies (> 
0.8 GeV) where fluxes are very similar one to each other. 

A cross comparison between Figure 15.41 and Figure 15.51 gives information on how cross 
section parameterizations (Kamae's and Tan's ones) changes the estimation of positron 
flux. As we observe, the differences among those are almost negligible if we compare them 
with the effect of propagation uncertainties. 

Uncertainties from nuclear cross section and CR ffux parameterizations show to be not 
so important as the uncertainties in the propagation. In any case, the former one would be 
important in an analysis which involves multimessengers correlations. 



5.5 Positron fraction analysis. 

The observable, that actual experiments are able to release flrst, is the positron fraction 
(Equation 14. 13|l in which the systematic errors are easier to manage and require less data 
processing than the flux determination case. The most controversial issue is the presence 
of an excess in the positron fraction, that was noted by Moskalenko et al. |T9 , when 
they derived predictions for the electron and positron fluxes. This excess has tried to be 
explained by many ways, for example, invoking contribution from DM annihilation (e.g. 
[iO] ) or contribution from nearby pulsars (e.g. [1121 [^ ). 
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In the later October 2008, the PAMELA experiment released the data for the positron 
fraction [H] , which would confirm a similar feature that HEAT experiment [SU] had slightly 
observed. Nevertheless, it is still unclear if this feature comes from a undiscovered source 
of positrons or maybe lay in other type of explanation. Another perspective is to think in 
the possibility that this feature is connected somehow to the electron component instead 
of the positron one. 

In the previous sections, we discussed and explained how to calculate the secondary 
positron flux. As well, we showed that current flux measurements are reproduced for same 
setting of the TZPM that reproduce B/C [H], proton and antiproton CR measurements [77] . 

In the work of Casadei et al. [53], they estimated the interstellar electron flux by com- 
bining the available data and producing a global fit. They performed on every dataset a 
demodulation, by inverting the modulation process described by Perko et al. [HI]. As well, 
to be able to demodulate the signal, they studied correlations between solar activity, other 
species of CR and information from neutron monitors [84] , in order to estimate the modu- 
lation parameter. The interstellar electron flux was modeled by a power-law function valid 
from 3 GeV up to 2 TeV: 



where 7 = -3.44 ±0.05 and N = (412.3 ± 22.8) m'^ sr'^ sec"! GeV"!. 

Figure 15 . 61 shows the positron fraction obtained with the two extreme cases for the 
electron flux, soft and hard. The soft electron flux has a power index of -3.54 (left) and the 
hard electron flux has one of -3.34 (right). These two situations correspond to a variation 
of 3(7 with respect to Casadei et al. reported values. 

In the case of a soft electron flux, a steeper behavior on the positron fraction is obtained. 
At high energies, the positron fraction tends to increase its value, and this is produced by 
the faster reduction in the electron flux. When we include the uncertainty band related 
to the propagation model, most of the experimental data goes inside this. On the other 
hand, in the case of a hard electron flux, the positron fraction in the high energy range 
gets reduced. This feature associated to a possible positron excess is supported better 
by the observations in the case of hard electron fluxes than in the case of soft electron fluxes. 

These extreme behaviors remark the important role of electrons in the positron fraction. 
As well, with the recent release of PAMELA data [H] the feature - above 10 GeV ~ has 
been observed, but the results at lower energies need to be analyzed carefully. 




(5.13) 
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Figure 5.6: Positron fraction as a function of the positron energy, for a soft and hard electron 
spectrum. Data are taken from CAPRICE [ST], HEAT [30], AMS 07], MASS ^ 
and PAMELA [41]. 



Experiment 


no [m"2 SJ--1 g-i GeV^i] 


7 


[GeV] 


xVd.o.f. 


AMS 


0.238 


3.86 


0.574 


0.152 


HEAT 


0.354 


4.03 


1.019 


0.237 


BETS 


0.2029 


3.201 


0.1122 


0.7985 



Table 5.3: Parameters values for best fits for electron flux of AMS [37] , HEAT [30] and BETS [33] 
experiments. 



An alternative analysis. 

A different approach is to analyze separately the electron fluxes of some experiments. We 
focused on the AMS [37], HEAT [30] and BETS [33] data, because those show different 
characteristics. For example, AMS and HEAT are more precise in the range of low-energy 
(< 30 GeV) instead of BETS which has measured the high-energy range giving valuable 
information about the IS flux. 



For these three experiments, we performed fits based on the Top-of- Atmosphere (TOA) 
electron data, using the functional form: 

1^10 GeV J U + log(10 GeV/0) ' > 

which depends on three parameters and it is clearly inspired by a modulated power-law 
function. The values of the parameters are given in Table 15.31 Notice that the value of 
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Figure 5.7: Demodulated (interstellar) electron flux versus electron energy. The flts of data 



from AMS [H7], HEAT [30] and BETS [31], have been demodulated according to 
modulation procedure used by Casadei et al. ^3]. The global fit obtained by Casadei 
et al. and the Strong et al. [39] flux estimation are shown as references. The AMS 
and HEAT fits are the ones with a behavior similar to Casadei et. al. On the other 
hand, BETS fit produce larger differences. 



X^/d.o.f. for each fit is lower that one, which indicates that the data can be nicely fitted 
with this functional form. 

Probably the logarithmic term cannot be explained by some physical phenomenon. How- 
ever, functions with that form behave like power-laws with variable power index. To clarify 
this point, the power index can be extracted from a function f{x) = a;'' doing: 



Applying the same method to x'^ \og{x/4)) , we obtain: 

power indexfx) = 7 + 1 — -, — i—r , (5.16) 
log(x/(/)) 

which shows that the power index changes slowly when x ^ (f> and it asymptotically 
converges to 7. 

Figure [577] shows the fits to electron data of AMS, HEAT, BETS and the Casadei et 
al. fit. In order to compare them, each flux was demodulated according to modulation 
parameters suggested by Casadei et al. (Table 13.21) . We notice that all curves exhibit 
similar behaviors except the one associated to BETS experiment, which present a harder 



power index(a;) 




power index ( a;) — 7 . 



(5.15) 
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spectra compared to the other experiments. Let us note that these smah differences may 
produce big effects into the positron fraction. 

We compute the positron fraction using the fitted fiuxes. Nevertheiess, we are not 
interested in processing the fluxes as Casadei et ai. did. Instead of that, just the positron 
flux is modulated, because the electron fluxes represent the TOA data which already 
contain modulation. 

Figure 15.81 shows the fits for each experiment and the correspondent positron fraction. 
In the plots for electron flux - left column - the la uncertainty bands (lighterblue) in the 
parameter determination are included. We can see that the fit for HEAT presents a narrow 
uncertainty band because its data are more accurate. The fit for BETS presents wider 
uncertainties at low-energy because there are not enough data in that range. The same 
situation occurs for the fit of AMS data, but in the high energy case. 

The right column of Figure [5~8l shows the positron fraction calculated with the positron 
flux obtained with the MED set (red line) and with the uncertainties in the propagation 
(darker blue band). Also, the flts for the data of AMS, HEAT and BETS, are used to 
calculate the fraction. The uncertainty band related to the parameter determination 
(lighter blue band) is included respect to the MED set situation. 

Depending on the experiment used to determinate the electron flux, we notice that the 
high energy range tail of the positron fraction may sizeably change. In the case of HEAT 
and using a modulated positron flux at 800 MV, the calculated positron fraction above 
10 GeV seems to be not as high as the experimental data, which means that a positron 
excess would be present. 

The same situation occurs when the AMS electron flux flt and the positron flux 
modulated at 600 MV are used to calculate the positron fraction. In this case, uncertainties 
in the parameter determination of the electron flux makes hard to determine if the excess 
is present or not. On the other hand, the excess is favored when the BETS electron fit is 
considered. 

The positron excess problem, that was noted as a deviation from Strong and Moskalenko 
predictions [39] , has attracted the attention of the scientific community. Most of the pos- 
sible solutions invoke the existence of sources like dark matter annihilation, primordial 
black holes and nearby pulsars. Nevertheless, the positron excess might not be exclusively 
related to an unknown contribution to the positron signal. This can be related to the 
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Figure 5.8: Electron flux s^^^+ (left) and positron fraction (right) versus energy. Each of the 
positron fraction has been calculated by using the specific data fit of HEAT [30] . 
AMS [37| and BETS [33], which are shown in the left column. The positron fractions 
were calculated using the positron flux for the MED set - red line - (See Table [472)) . 
The uncertainties related to the propagation (darker blue) and the la uncertainty 
(lighter blue) band for each flt are included. Experimental data were taken from 
CAPRICE tSlj, HEAT AMS t37j, MASS 02], PAMELA [H] and ATIC-2 [35] 
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electron signal as well. We showed how different data sets from some single experiments 
may sizeably modify the predictions for the positron fraction. This docs not discard the 
contribution of new sources: on the contrary, this motivates to study deeper and consistent 
the full picture, including the theoretical estimation of the electron flux. 
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Chapter 6 



Conclusions 



In this thesis, we studied different aspects of positron cosmic-ray physics. The main 
aims were the study of signals related to galactic dark matter annihilation and secondary 
positrons. The propagation of positrons were studied, as well as the effect of uncertainties 
related to the cosmic-ray propagation. We described and estimated the impact of uncer- 
tainties on the positron flux. 

In the galactic enviroment, the positron production takes place in many different ways. 
We described and studied two of the most important processes (Chapter [2|). The first one 
is the production of positrons and electrons in nuclear spallations in which positrons and 
electrons are mainly produced through the decay chains of charged pious and kaons. The 
second mechanism of production is related to dark matter annihiliation. We studied and 
characterized positron and electron energy spectra obtained from neutral initial states. 
The method consists in the production of particlc-antiparticle pairs in which the energy 
spectra are generated from the simulation of decay chains and hadronization process 
based on Lund's PYTHIA Montecarlo For the case of more complex intital states, we 
calculated them in terms of former spectra using the rules from the Standard Model and 
the Two Higgs Doublet Model. 

The propagation of positrons, and in general of any type of cosmic-rays, was treated 
according to the Two-Zone Propagation Model [Bl]. We solved the positron-electron 
transport equation using the Green fvmction method and we defined the concept of 
Halo-Function, which is the response of the propagation domain to the distribution of 
sources. We calculated the halo functions for three dark matter distributions - cored 
isothermal [50] j NFW [HI] and Moore's profiles - and for the source distribution of 
secondary positron production. The halo-function concept naturally includes a definition 
of the diffusion lenght Ad (Equation I3.28|l . which relates the propagation distances with 
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energy evolution of cosmic-rays. When energy losses related to inverse Compton scattering 
and synchroton radiation dominate the transport equation, the halo-function depends 
exclusively of Ad and L. This fact improves the speed in the calculations of the positron 
signal, as well, it makes accesible to realize scans over the propagation parameters. 

The geometry and boundaries of the propagation zone are important and affect 
the solutions of the transport equation. The Two-Zone Propagation Model considers a 
propagation zone described by a cylinder with radius 20 kpc and thickness 2L, which 
describes the extension of galactic magnetic fields reponsables of cosmic-ray diffusion. 
We calculated Green functions for cartesian and cylindrical coordinates with and without 
boundary conditions. We found that the presence of boundaries speeds up the leaking of 
particles from the propagation zone, that is reflected in the halo function vanishing when 
Ad is larger than the distance to the closest boundary. 

The propagation uncertainties on the halo function have been calculated for the 
1,600 different cosmic ray propagation models that have been found compatible [64] with 
the B/C measurements. These uncertainties are strongly dependent on the source and 
detection energies, Ss and e. As gets close to e, we observe that each uncertainty domain 
shrinks. In that regime, the diffusion length Ad is very small and the positron horizon 
probes only the solar neighborhood. In the opposite case, the uncertainty can be as large 
as one order of magnitude or even more. As positrons originate further from the Earth, 
the details of galactic propagation become more important in the determination of the 
positron flux. On the contrary, high-energy positrons are produced locally and the halo 
function / becomes unity whatever the astrophysical parameters. 

In the case of dark matter annihilation, we inspected specifically the positron fluxes 
for a 100 GeV DM particle annihilating into four typical channels: direct production of 
positrons, r+r^, W^W~ and bb pair. Uncertainties due to propagation on the positron 
flux are one order of magnitude at 1 GeV and a factor of two at 10 GeV and above. We find 
an increasing uncertainty for harder source spectra, heavier DM and steeper proflles. 

The comparison with current data shows that the possible HEAT excess is reproduced 
for DM annihilating mostly into gauge bosons or directly into a positron-electron pair, 
and the agreement is not limited by the astrophysical uncertainties. A boost factor of 10 
is enough to obtain an excellent agreement between the measurements and the median 
flux, for a 100 GeV DM particle. A smaller value would be required for a flux at the upper 
envelope of the uncertainty band. 
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We have finally drawn prospects for two interesting 3-year flight space missions, like 
PAMELA, already in operation and which has recently released its first data [H], and the 
future AMS-02. We reach the remarkable conclusion that not only will PAMELA have the 
capability to disentangle the signal from the background, but it will also distinguish among 
different astrophysical models, especially for hard spectra. For AMS-02 the possibility to 
disentangle the signal from the background is also clearly manifest. We also wish to remind 
that improved experimental results on cosmic ray nuclei, expecially on the B/C ratio, 
will be instrumental to improve the determination of the parameters of the propagation 
models, and will therefore lead to sharper theoretical predictions. This in turn will lead to 
a more refined comparison with the experimental data on the positron flux. Moreover, a 
good determination of the unstable/stable nuclei abundances like the ^°Be/^Be ratio could 
shed some light on the local environment, which is certainly mostly relevant to the positrons. 

Furthermore, we studied the secondary positron signal and we have compared the vari- 
ous models available for the interstellar secondary positron production. It has been shown 
that more positrons are expected when the proton flux from Shikaze et al. [47] is used, as 
compared to the case proposed by Donato et al. [S5] . Moreover, for a given proton flux, the 
three positron production cross sections we have considered produce different result: below 
a few GeV, the parameterization of Badhwar et al. [48j gives more positrons, whereas above 
a few GeV, the model of Tan and Ng |50] predicts a larger positron production. At any 
energy, the parameterization of Kamae et al. [56^ produces the smallest amount of positrons. 

We showed that varying the diffusion parameters has not the same effect as for primary 
positrons [57]. For exotic positrons created in the Dark Matter halo, the thickness of the 
slab 2L was the most relevant parameter because the increase of the diffusion zone implies 
the increase of the number of sources, whereas for secondary positrons - which are created 
in the galactic disk only - the most relevant parameter is the diffusion constant i^o- 
Therefore we expect the sets of parameters that basically maximize the primary positron 
flux to minimize the flux of secondary positrons, and vice versa. 

We also showed that, because of energy losses during propagation, most of the positrons 
detected at the Earth have been created in the nearby 2 kpc: this is the reason why we 
could safely neglect the variation of proton flux in the galaxy. 

Finally, and this is our most important result, our estimation of the positron flux is 
compatible with all available data. This does not mean that there is no exotic positron 
contribution, since we have not tried to fit the data with a single diffusion model. However 
this shows that one should be cautious before claiming that there is any excess in present 
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data. As regards a possible excess in the positron fraction, we have also clearly shown 
that the electron flux actually plays a role which is as important as that of positrons. 
This might sound tautological because the positron fraction is no more than a ratio, but 
so many energy has been involved in this issue that we again strongly stress this simple fact. 

The just released PAMELA data [H] show a clear rise of the positron fraction for 
energies above 10 GeV. From our analysis, whose objective is an accurate determination of 
the positron flux, we derive the conclusion that the presence of an excess is clear in the 
case of a hard electron spectrum, while for a soft electron spectrum the rise of the positron 
fraction may be explained by the standard secondary production. Scanning over the various 
parameters at stake, we in fact find in general the PAMELA measurements in excess of 
what a pure secondary component would yield. Nevertheless, if the electron spectrum 
is soft most of the PAMELA data points are aligned with our prediction [105] . Notice 
that, also in that case, the two last energy bins feature an increase, but the experimental 
uncertainties are large there and a presence of an excess is, in this case, currently not 
statisticUy significant. 

In this thesis, we have thus presented the methods and the practical tools to evaluate 
the primary positron fluxes in detailed propagation models. We have provided careful 
estimations of the underlying uncertainties and shown the extraordinary potentials of 
already running, or near to come, space detectors. 

More insight in these issues will therefore require, from the theoretical side, a revised 
understanding also of the electron flux, including the determination of its uncertainties, 
and from the experimental side, the separate provision of the electron and positron fluxes, 
in order to better compare theoretical predictions with the data. In addition, the upcoming 
data on cosmic rays above 10 GeV will allow us to reduce considerably the theoretical 
uncertainties on all cosmic ray fluxes and help us elucidate the experimental status of the 
so-called excess in the positron spectrum. 
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Appendix A 

The boosted spectra method 



Particle distribution functions are easier to calculate in the center of mass 
frame. However, this situation is not the most common one. Doing some 
kinematical analysis, the distribution can be rewritten for different reference 
frames in terms of new kinematical variables. 

The first aim of this method is to transform the energy (spatial) distribution function 
from a specific reference frame into another. In Chapter [H it is used for calculating the 
positron (electron) energy spectra in the case of dark matter annihilation processes. 

Let us start taking an angular and energy dependent spectrum, 

/(e,cos0) = — . (A.l) 

Under change of frames, the spectrum shape would change although it keeps constant its 
integral. 

The energy Lorentz transformation between two frames is: 

e' = 7(e-/3pcos6') , (A.2) 

which gives restrictions on how energies will redistribute. All physical configuration, in the 
prime frame, are considered when the spectrum is integrated angular and in energy, 

/'(e') = \ de d{cose) /(e,cos6i) (5(r(e';£,cos6i)) , (A.3) 



where r is a restriction obtained from Lorentz transformations. Angular integration gives: 

= ^ /-/(£, cosr), (A.4) 
IP J P 
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Figure A.l: Schematic graphics E' vs E. It shows how to determine integration bounds 

for a fixed energy e' (in green). Maximum and minimum energies for the prime frame 
in terms of original frame energy (in red) are also shown. 



where, 



cosr 



7£ — e 
7/3 p 



(A.5) 



is the angle respect to the boost direction, which links energies between frames. A simple 
case is an isotropic spectrum, where the boosted spectrum is: 

de 



fLie') - ^ y 7 fUs) 



(A.6) 



A very important point is about the integration limits on the original frame energy. 
Those limits come from inverse Lorcntz transformations and from the maximum energy 
(^max) allowed in the original frame (Figure [A. Then the range on energy in the original 
frame, which contributes to the spectrum at a fixed energy e' is: 

emax = min (^Smax, 7(e' + /? P')) (A.7) 
Emin = max (^m, j{e' - /3 p')) 

The maximum energy in the prime frame, 

Smax = 7(emax + P Pmax) , (A.8) 
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is obtained by transforming Emax- That provides a complete description for transforming 
any kind of spectrum. 
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Polarized Muon decay 



The Standard Model and Fermi theory give a precise description of muon 
decays. A common situation considers unpolarized muons decaying into 
electrons and neutrinos. However, there are processes where muons are 
produced fully polarized, as a consequence of left-right asymmetry present in 
electroweak interactions. 

A production mechanism of electrons and positron is possible through the muon decay. 
Muons decay emitting two neutrinos and a positron or an electron (depending of muon's 
electrical charge). Usually muons are produced unpolarized, however when those come 
from a spin-0 state, they become fully polarized. That process modify the final electron or 
positron energy spectra. The polarized muon state induce an angular distribution respect 
to muon's propagation direction. 

Positrons and electron produced in nuclear spallations are more affected by the muon's 
polarization efi'ect, because the amount of charged pions and charged kaons produced 
enhance this effect and produce deviations from the unpolarized case. 

In the muon rest frame, the electron energy and angular spectrum is [2]: 



where x = 2e/m^ and y — cos 6. Those are related with the electron energy and angle 
respect to the muon polarization vector. 

It is important to remark that this spectrum is normalized to unity. 



dxdy 



fix, y) =x^{3-2x± y{2x - 1)) - fi{x, y) ± /^(a;, y) , 



(B.l) 




(B.2) 
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in other words, it can be directly related to a probability density. 

A property present in the muon rest frame is that the anisotropic term cancels when 
angular integration is performed, 

/(x) - J ^dy fix, y) = 2x^ (3 - 2x) . (B.3) 

When muons are in motion, angular distribution becomes important and we need to con- 
sider a general approach. To solve that we used the boosted spectrum formulae (Appendix [1]). 

In this context, the boosted energy spectrum is: 
where 

is related with physical angle which gives the correct energy x' . 

The integration limits are obtained from inverse Lorentz transformation: 

f x'/x'^^ 0<x'<x^i, 

X I 

where 

x'=l{^ + P) and = 7(1 - /3) ■ (B.7) 



For simplicity, electron boosted spectrum can be divided into an isotropic and an 
anisotropic term. 



f{x')^fAx')±fAx'). (B.8 



After the integration on x, the isotropic part becomes: 

rj{x')^lx'HQ-2x'j{3 + (3')] , (B.9) 
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when < x' < x^j^, and 



when <x' < x'^^^. 



Also the anisotropic part is sphtted into two expressions: 

2 



f'A{x') = lx'^imx'^-i) (B.ll) 



and 



f, _ l + x"x'i[2x'x',^{l + W)-^l + W)] ,^ , 

for same the ranges than the isotropic case. 

Let us clarify that we work in terms of dimensionless quantities. In order to get the 
energy spectrum in proper units, we need to do: 

which has units of [e~^]. 



B.l Accumulated probability 

The accumulated probability is a useful tool to be used in Montecarlo programs because it 
can be used for event generation. By definition, 

F{x') = r f{t')dt' , (B.14) 
Jo 

gives a monotonic function with values which go from to 1 . Computing the accumulated 
probability for the isotropic and the anisotropic terms, we obtain that: 

Fi{x') = la;'S[6-x'(3 + /32)7], (B.15) 

Fa{x') = ^x'V[2A-l], (B.16) 
for < x' < x'^f^. And with values at x'^^i- 

3 + /3[3 + /?(/3-7)] 2/3(1-/3)^ 
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Distribution function for x' between x'^^ and Xg^p are: 

Fi{x') = i^/«h) + ~ ^ch) + <h[<h(^^'^- <t) - 3(a;'^ - x'l)] ^^^^^^ 



To produce electron (positron) energy spectrum in a Montecarlo, it is necessary to invert 
these function. This can be obtained by using numerical invertion methods. 
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Electrons and Positrons energy 
spectra from Pion decay 



The principal source of fully polarized muons are the charged pions. Those 
are spin-0 particles which are abundantly produced in nuclear scattering 
processes. 

Charged pions decay mainly into muons with a probability of 99.98770 ± 0.00004% [5]. 
Nevertheless, muons are produced fully polarized due to pions are spin-0 particles (See 
Figure \CA^ . 

Depending on the pion charge, electrons or positrons are more or less energetic with 
respect to unpolarized situation. This phenomenon should be taken into account in order 
to understand possible asymmetries present in this type of CR. 

In the pion's rest frame, all muons come from a two-body decay mode. In this frame, 



Figure C.l: Schematic plot to illustrate how the charged pions decay produces polarized muons. 

Blue arrow represent left-handed polarization and red arrows right-handed polarized 
particles 
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they are monoenergetic, and then, the energy spectrum associated is: 

^ = ^(eA. - £*) = //.,7r(e/.,TO^) , (C.l) 

where e* is the muon energy in pion rest frame and its value is obtained from energy- 
momenta conservation, 

ml — mt 



The electron spectrum can be obtained as a convolution between pion decay - into 
muons iffi^Ti) - and muon decay - into electron (fe.fi), 

ye,7r (^e ; ^7r) — ^^fi //^,7r (^/i ; ^tt) fe,^(^^e:^fi) • (^•'^) 



Muons are monoenergetically produced when pions are at rest f Equation IC.2|) . In this 
limit the electron spectrum is simply: 

fe,^{£e,m„) = fe,i,{ee,£l) . (C.4) 

When pions arc in motion, the boosted spectra method are applied to the muon spectrum 
. In this situation, the boost factor depends on pion energy 7^ — £^/m,r, and then, the 
muon spectrum is: 

\ ^min ^ ^ ^ ^max 

ff,,^{£f,, £7r) ^ < , (C.5) 



other cases 

where energy bounds are got from boosting e*, 

e'r^-f^i^* +(^-P*) and £-- = 7,(£*-/3,p*) . (C.6) 
Finally, the electron spectrum is calculated through convolution of partial spectra: 

/e,7r (^ej £ir) = T: 3 - / ds fe,fi{£e,£fj.) , (C-7) 

where the information on muon's polarization is inside /e,/x- 
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The maximum electron energy for a pion decay can be calculated from composition of 
Lorentz transformations from muon and pion rest frames, 

777 

= 7.(1 + 7r'(l + f3rn^ ■ (C-8) 



c.l Analytic formulae for polarized muons 

The electron spectrum from pion decay depends on the integral on muon energy. As it was 
seen in AppendixEl electron spectrum from polarized muons can be described by composition 
of isotropic and anisotropic terms, 

feA^e,e^) = — {f'j{x')±fAx')) , (C.9) 

where x' — 2ee/m^ and prime denotes that those were calculated for a moving muon with 
boost factor of 7 = e^/tti^. 

The integration on muon energy can be transformed in terms of the muon's 7 factor, 

de^ feA^e.e^) = 2 T'™ rf7 {f'j{x') ± f^x')) , (C.IO) 

and then work in terms of dimcnsionlcss quantities. 

An important point is about different behaviors present on /'. To work in a abstract 
way, it has been renamed the formulae for muon decay (showed in Appendix |b|) for each 
range of x' , 

0<a;'<<h^/{(x') ; < ^' < <up ^ /al^') ■ (Cll) 



When 7 varies, x'^-^ and xi,^p vary too. Then 

7- = (C.12) 

marks the critical point where /' changes behavior. Depending on x' value, the integral 
just considers one or both functions (Figure [C.2|) . 

For x' <1 and 7„iin < 7* < 7max: 

d7/'= r d^f[ + r"^'^d-,f!, (C.13) 
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Litieros' PhD thesis f 




Figure C.2: Schematic plot x' versus 7. Integral limits on muon's 7 and frontiers between zones 
(red line), where spectrum changes behavior, are shown. 

Otherwise, it is the integral which corresponds to the proper sector. 
For x' > 1: 

(C.14) 



J max(7* ,7min) 

where there is contribution from because region below 7* are not allowed. 



Performing the integral on the muon's 7 factor over isotropic and anisotropic terms, 
a set of indefinite integrals are obtained. Where they should be evaluated in the proper 
integration limits. 



The integration of isotropic term gives, 

„2r 



Fnixn) = ^(277 + 4a;(3-472)) (C.15) 
Fi2{x,^) = ^(x2{27x^h7-4x7(/?-3 + 44h7)} + 151og(x^„p)) , (C.16) 
and the anisotropic term is: 

Fai{x,i) = y(7'/3(16a;7/3'-9) + 91og(x^„p)) (C.17) 

,(x,7) = ^(2x2{9x;i^7-8a;7(2a;;h + 2/3-3) + 91ogK,p)} (C.18) 

+ (3 - 9a;2) \og{l3'' -f'^) - Gx^ log(lii)) 

7 — 1 / 
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Finally, we have obtained all the ingredients to compute a realistic electron spectrum. 
Those are been used to replace muon decay subroutines in the Lund's PYTHIA montecarlo 
generator [61]. 
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Appendix D 



Electrons and Positrons energy 
spectra from Kaon decay 



The second in importance source of fully polarized muons are the charged 
kaons. Those are spin-0 particles which are abundantly produced in nuclear 
scattering processes, as charged pions. 

The electron spectrum from charged Kaon decay can be more difficult to manage due 
to the various decay modes of Kaons (Table ID. l|l . These modes are dominated by two body 
modes with a probability of ~ 85%. The remained modes are mainly composed by three 
body ones. As in pion decay, muons are produced fully polarized. 

To calculate a complete electron spectrum, each decay mode can be treated separately 
and weighted according to its own probability. 



In this way, each mode is independently studied. 



D.l Kaon two body modes 

Two body decay modes are easily described by analytical expressions. The electron 
(positron) spectrum, for isotropic two-body modes, is: 




(D.l) 



mode 




(D.2) 
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K decay mode 


Probability % 


(Tie) produced 




63.44 ±0.14 


1 




20.92 ± 0.12 


1 




5.590 ±0.031 


2 




4.98 ±0.07 


1 


tt" /i 


3.32 ±0.06 


1 




1.757 ±0.024 


1 


Total 


100.0 ±0.445 


1.06 



Table D . 1 : Most important charge kaon modes, probability and number of electron associated 
to each of these [2]. 



which is for Kaon at rest. On the other hand, for Kaons in motion with boost factor 
7 = Ek/sk, the spectrum is: 

fe,K{£e,£K) = „ „ ^ / dsx fe,x{£e,£x) , (D.3) 

where particle X refers muons or charged pions or any two-body mode. The integration 
limits are calculated from e'^ boosted in the backward and forward direction respect to the 
boost direction, 

= 7(ex + Px) , eT' = l{ex - P Px) , (D-4) 
Where p* depends on which two-body modes is involved. For the mode K~ — > iiP^, it is: 

(D.5) 



and for — > tt^tt^ , 



" 2^ V + " "''^"''^ ^ ml^m]^ . (D.6) 

Note that when Kaons decays into Pions, those modes are more complex than modes with 
muons. Due to that, it is needed to perform numerical integration. 



D.2 Kaon three body decay 

The next important contribution to electron spectrum are composed by three body modes. 
124 



D.2. KAON THREE BODY DECAY 



D.2.1 Mode 7r° e 



This is a mode where energy spectra for positron and electron are equal. For an isotropic 
production, energy spectrum at Kaon rest frame is: 



where Uq = 54.14 GeV-^, ai = 0.2468 GeV and E^,^^ = 0.2283 GeV. 



(D.7) 



The boosted version of this spectrum is: 



/ttci^ (^e? ^/c) — 2'y/? ^ 



where 7 = EK/friK and 



7(1 + 

Be. 



7(1-/3) 
E 



aJmin + (ai - -Bmax) log( ^°'''''_ ) ) , 



< ee < 7(1 - fi)Era^^ 

7(1 - /3)^„,ax < £e < 7(1 + /3)^n 



(D.8) 

(D.9) 
(D.IO) 



D.2.2 7r° /X 

For this mode, it is needed to calculate the muon energy spectrum, where the maximum 
muon energy at kaon rest frame is: 



and the spectrum: 



s\ + ml - mlo 
<,ma. = = 0.2395 GeV, 



e / \ A (^/i,max ^m) 

htii^,K{eti,eK) = A p^- — — . 

oo — ai£^ — 



(D.ll) 



(D.12) 



A = 1374.86 GeV^i (D.13) 
a„ = 0.5046 GeV^ (D.14) 
ai = 1 GeV (D.15) 



The boosted energy spectrum is: 



fniJ.i^,K{Sn,S:K) = {S'7r/ii/(emax) — ^^(^^(einin)} (D.16) 
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where QT^^uit) is an auxiliary function defined: 



where 



Ao = V4ao + a? (D.18) 



In order to get the electron spectrum, it is needed to convolute this spectrum with the 
electron spectrum due to muon decay. 

max 

fe,K{£e,£K) = / de ^ Itt iiu,k{£ £ k) /e,^i (^e , , (D.19) 

In this case: 

£r'=^«max + /3p;,max) ; sf^'^m^ (D.20) 



D.2.3 11-71-11+ 

This mode produces a couple of electrons and one positron and muon polarization effect 
are also present. As it is known, electron spectrum can be got from convolution between 
the pion spectrum from kaon decay and electron spectrum from pion decay (Appendix [c|) , 

/e,_ff(ee, ~^ J '^^'^ f3n,K{£ir,£K) /e,7r (^e , Ett) , (D.21) 

and it is convenient to work just on the 3 pions energy spectrum. Note that the factor 2 
is related with the electron multiplicity. 

At kaon rest frame, the maximum pion energy is: 

<nax - = 0.18758 GeV (D.22) 

2eK 

and the pion energy spectrum can be written as: 

f37T.K{£yr,£K) = ^ Ptt (Ett " m^f (Sniax " Syr)^^^ (D.23) 

where the parameters are: 

A = 2002.19 GeV"2-5 ; S = 0.05 (D.24) 
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for the actual measured masses. 

The boosted version is obtained through: 

/37r,K(£7r,£x) = / ~ f3-n,K{e,mK) , (D.25) 

T/^ J €min 

where 

emin = liSn-PPTr) (D.26) 

l{e„+PPn)) (D.27) 

and then, it can be written in a analytical form: 

f3Tr,K{sTT,eK) — {S'37r(emax) — fi'37r(emin) } (D.28) 

by using the auxiliary function gsj^' 

53.(e)=t(e)v/7^f;^V 2fJs + 1;S-^;S + 2;^) (D.29) 

where 2-Fi(a; b; c; z) is a hypergeometric function, i(e) and Tmax are the kinetic energy and 
the maximum kinetic energy for charged pions. 



D.2.4 7r-7r°7r° 

This mode is very similar to the previous seen. Kinematics restrictions are slightly different. 

£max = — = 0-1926 GeV (D.30) 

The same parameterization is good enough to describe the energy spectrum at kaon 
rest frame: 



where the parameters are: 

A = 1624.47 GeV~2-^ ; S = 0.065 (D.32) 

And the boosted version is: 

fw±2-n>>,K{£w±,^K) = {537r(emax) - 537r (emin) } (D.33) 
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where 



Cmin = ^iSn-Pp^) (D.34) 

Emax = min ( e^^^, j{e„ + /? Ptt) ) (D.35) 



and (?37r(e) is the same auxihary function used for the mode with three charged pions. 
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Method of Averages: A method 
to produce accurate histograms 

The Method of Averages (MoA) is an algorithm which has been developed 
to estimate the proper binning rule to create a histogram. This method uses 
directly the data as an input and produces a binning rule where narrow bins 
are placed in high-density zones and wider bins in less dense zones. 

There are many criteria to estimate the number of bins necessary to make a histogram 
fTable |E~T|1 . Nevertheless, in most of cases, there are big assumptions about how the data 
are distributed. As an example, Scott's choice [114] makes the assumption in which data 
follow a gaussian distribution, where bin widths are proportional to the standard deviation 
of the sample. 



Model 


number of bins 


Sturges' formula |113j 
Scott's choice [114] 
Freedman-Diaconis' choice [115] 


log2{n) + 1 

-^max -^min l 

n ' 

3.5 a 

3^max ^min i /3 
ft 

2 IQR(a;) 



Table E. 1 : Number of bins estimation formulas for Sturges, Scott and Freedman-Diaconis models 
|113l 11141 1115j . For all models, n is the number of data, x^ax / min are extreme values 
present, a is the standard deviation and IQR(a;) is the interquartile range function. 
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DATA 




c> x„,„') pass-1 






C (•'^> ) 


( {x) ) pass-2 


(W) CW) 

mm 


( jx) ) ( (x) ) pass-3 

414 







Figure E . 1 : Schematic view of creation of MoA-rule and how data splitting is performed. 

E.l Description 

The key in Mo A is to analyze the sample in a simple way. For that, we classify extracted 
quantities in term of number of passes, or times in which data were used. After reading the 
data once, the pass-1 quantities are the extreme values, a;i„ax/minj and the average {x)pi. 
The pass-1 average gives information about the tendency point; for example, if the average 
is close to Xinin, that means data are roughly concentrated around the minimum. However, 
if the pass-1 average is close to the mid point between Xmax and Xnim, we expect that 
data are spread in the whole range. Of course, the real information is inside the standard 
deviation, but that belongs to the pass-2 quantities group. 

To get an accurate binning rule, we need to read data many times. To reduce the 
screening effect due to dense-data zones, before going to a next pass, we split the sample 
into subsets. Those are delimited by averages obtained in previous passes (Figure IE. l|) . For 
example, to compute pass-2 averages, we need to compute the averages of data in ranges 

(•^min : 

(x)pi) and ((x)pi,Xmax)- 

At the end of p passes, we obtain a rule composed by umoA — 2^ — 1 bins. The choice 
of p depends on the desiderate number of bins (nbin)- It is recommendable to use a value 
of p that satisfies: 

nMoA < , (E.l) 
to reduce the risk to skip one of the averages, which are very important point in our method. 

Starting from the rule obtained with MoA, we generate the final rule using linear in- 
terpolation among MoA-rule points (Figure lE.2p . The interpolation is done assuming that 
MoA-rule points are homogeneously and monotonically distributed in the range (0, 1) of a 
parameter i, where t = corresponds to Xmin and i = 1 to Xmax- For example, a MoA-rule 
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Lineros' PhD thesis (2008) 




Figure E . 2 : Example of rule created with MoA and the composition of the final rule by using lin- 
ear interpolation. In this figure: red dots are the averages and red crosses correspond 
to extreme values. Note that final rule presents a non homogeneous distribution. 

composed by n -I- 1 elements is mapped into n + I values of t according to: 



„MoA 



.MoA 



- y i eZ / <i <n 
n 



where the extremes are Xq 



MoA 



and Xq 



MoA 



(E.2) 



To get the final rule with ribin bins, a homogeneous discretization on t is performed and 
then the linear interpolation between MoA-rule points is easily managed. 



The MoA can be easily tested and implemented. For example in Figure IE. 31 the distri- 
bution for generating random events is: 

p{x) = ax~^ +I3exp - Wf^ +(3exp - lOOf^ , (E.3) 

which is a composition of a power-law distribution and two gaussian distributions, one 
centered in a; = 10 and the other in a; = 100 with the same standard deviation of 1, finally 
a and P are just normalization constants. Usually, power-law distributions are better 
described when a uniform-log rule is used. On the other hand, gaussian distributions are 
usually well described with an uniform-linear rule. When both distribution are mixed, 
there is not an efficient binning rule based on generic rules. After both use the MoA to 
analyze the generated data, the binning rule produced describes well enough Gaussian and 
the power-law distributions. 
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X 



Figure E . 3 : Histograms made with a sample of 10 events and 60 bins. An analytical distribution, 
MoA and uniform log-bin rule are plotted. The analytical distribution is composed 
by three known distributions, a power-law (oc s"^) with two gaussian located in 10 
and 100, all distribution have same weight. Binning rule obtained from MoA produce 
a more accurate solution than the uniform log rule, especially in the regions close to 
gaussians. 



E.2 Histogram creation 

Once the binning is produced, the next step is to scan the data and fill each bin of a 
histogram. To take advantage of the process of counting, we compute data-averages for 
each bin, 

/ \x- - f An\ "events e (a;j, Xi+i) 



This process increase the accuracy in the histogram because the tendency points are 
computed for each bin instead of the mid point of the bins, which is the standard way. 

This method gives good results in case of many combined distributions. Especially, this 
can be applied to produce particle energy spectra as in the case of particle decay chains, 
where due to particle boost energy spectra are deformed. 

The MoA was used in the production of histograms related to DM annihilation (See 
Chapter m . We generate the data with Lund's PYTHIA for several types of particles jets. 
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superimposed on each other |6C 

4.3 In each panel, the halo integral / is plotted as a function of the positron 
injection energy Eg whereas the energy e at the Earth is fixed at 10 GeV. The 
galactic DM halo profiles of Table 14. II are featured. The curves labeled as 
MED correspond to the choice of cosmic ray propagation parameters which 
best- fit the B/C ratio [SI]. The MAX and MIN configurations correspond 
to the cases which were identified to produce the maximal and minimal DM 
antiproton fluxes [77] , while the entire colored band corresponds to the com- 
plete set of propagation models compatible with the B/C analysis [53]. The 
MAX, MED and MIN configuration are given in Table [3TT1 |61 

4.4 Same plot as Figure [4~3] the Moore DM profile has been selected. Four values 
of the positron detection energy e have been assumed. The fiag-like struc- 
ture of this figure results from the widening of the uncertainty band as the 
detection energy e is decreased [62 

4.5 Diffusion length Ad versus injection energy. The line at 5 kpc has been traced 
to denote that Ad bigger than this limit represent the range where the central 
DM overdensities start to appear. The two extreme B/C analysis configura- 
tions and the best-fit MED (Table 13 . l|l have been used to calculated their 
associated Ad , showing how the MAX configuration covers longer distances 
instead of MED and MIN 63 
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4.6 Positron flux i?^$g+ versus the positron energy e, for a DM particle with a 
mass of 100 GeV and for a NFW profile (Table |47T|) . The four panels refer to 
different annihilation final states : direct e+e" production (top left), bb (top 
right), W^W^ (bottom left) and r+r^ (bottom right). In each panel, the 
thick solid [red] curve refers to the best-fit choice (MED) of the astrophysical 
parameters. The upper [blue] and lower [green] thin solid lines correspond 
respectively to the astrophysical configurations which provide here the maxi- 
mal (Ml) and minimal (M2) flux - though only for energies above a few GeV 
in the case of (Ml). The colored [yellow] area features the total uncertainty 
band arising from positron propagation 



4.8 Positron flux e^$e+ versus the positron energy e, for a DM particle mass of 
100 GeV and for different halo density profiles : cored isothermal sphere [90] 
(left panels), NFW [HI] (central panels) and Moore (right panels) 
(Table |4~T|) . The upper and lower rows correspond respectively to a bb and 
W^W~ annihilation channel. In each panel, the thick solid [red] curve refers 
to the best-fit choice (MED) of the astrophysical parameters. The upper 
[blue] and lower [green] thin solid lines stand for the astrophysical configura- 
tions Ml and M2 of Table 14.21 The colored [yellow] area indicates the total 
uncertainty band arising from positron propagation 

4.9 For fixed values of the detected energy e, the uncertainty bands on the 
positron flux e^m^$e+ are shown as a function of the mass of the DM 
particle. The energies considered in the flgure are e = 3, 10, 30 and 100 GeV. 
Each band refers to one of those values and starts at — e 



4.10 Positron fraction jieT + e"*") versus the positron detection energy e. No- 
tations are as in Figure 14.61 In each panel, the thin [brown] solid line stands 
for the background [301 ESj whereas the thick solid [red] curve refers to the 
total positron flux where the signal is calculated with the best-fit choice 
(MED) of the astrophysical parameters. Experimental data from HEAT 
AMS [361 [37], CAPRICE 01] and MASS [32] are also plotted 



4.11 Same plot as in Figure |4. 101 but with a mass of the DM particle of 500 GeV. 

4.12 Positron flux i?^$g+ (not fraction) versus the positron energy e, for a 500 GeV 
DM particle. Notations are the same as in Figure 14. 101 Experimental data 
from HEAT [SU], AMS [SHIST], CAPRICE [21] and MASS [32] are pfotted. 
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4.7 Same plot as in Figure |4~6] but with a DM particle mass of 500 GeV [67 
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4.13 Predictions for PAMELA for a 3-year mission. The positron fraction 
e+/(e~ + e+) and its statistical uncertainty are pfotted against the positron 
energy e for a 100 GeV DM particle and a NFW profile. Notations are the 
same as in Figure 14. 101 The thick solid curves refer respectively to the total 
positron flux calculated with the Ml (upper [blue]), MED (median [red]) and 

M2 (lower [green]) sets of propagation parameters |78 

4.14 Predictions for AMS-02 for a 3-year mission. The positron flux £^^e+ ^^nd 
its statistical uncertainty are featured as a function of the positron energy 
£ for a 500 GeV DM species and a NFW profile. Notations are the same as 
in Figure 14. 121 The thick solid curves refer respectively to the total positron 
flux calculated with the Ml (upper [blue]), MED (median [red]) and M2 
(lower [green] ) sets of propagation parameters [SC 

5.1 Positron source term e^'^^q versus positron energy |85 

5.2 Halo function for the secondary positron source term versus the diffusion 
length Ad (left). Corrective function /c versus Xu/L^- Both the numerically 
calculated function and the analytical fit are plotted (right) |86 

5.3 Interstellar electron and positron fluxes e^^e for the MED set (Table [3TT|) ver- 
sus energy. The curves correspond to fluxes calculated from Nuclei CR inter- 
actions with the ISM. Each curve represents a different nuclear cross section 
and Nuclei CR flux parameterization. Also the Strong et al. f Equation 14. 9|) 
flux parameterizations are shown. The uncertainty band related to those pa- 
rameterization is plotted (yellow band) as well |88 

5.4 Secondary positron flux £^'^^^+ versus positron energy. Positron fluxes were 
calculated using proton and alpha CR fluxes from Donato et al. [65] and 
Shikaze et al. [171 with the Kamae et al. cross section parameterization [SS]. |9C 

5.5 Secondary positron flux e'^'^$e+ versus positron energy. Positron fluxes were 
calculated using proton and alpha CR fluxes from Donato et al. [5S] and 
Shikaze et al. ^\ with the Tan et al. cross section parameterization [SU] ... |9C 

5.6 Positron fraction as a function of the positron energy, for a soft and hard elec- 
tron spectrum. Data are taken from CAPRICE [31 , HEAT [30 , AMS [37], 
MASS [32] and PAMELA [H] m 

5.7 Demodulated (interstellar) electron flux versus electron energy. The fits of 
data from AMS :37J, HEAT [30J and BETS [33], have been demodulated 
according to modulation procedure used by Casadei et al. [13] . The global fit 
obtained by Casadei et al. and the Strong et al. [32] flux estimation are shown 
as references. The AMS and HEAT fits are the ones with a behavior similar 

to Casadei et. al. On the other hand, BETS fit produce larger differences. . |94 
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5.8 Electron flux £^<i>e+ (left) and positron fraction (right) versus energy. Each 
of the positron fraction has been calculated by using the specific data fit of 
HEAT [30 , AMS [37] and BETS [33, which are shown in the left column. 
The positron fractions were calculated using the positron flux for the MED 
set - red line - (See Table |4T2|) . The uncertainties related to the propagation 
(darker blue) and the Icr uncertainty (lighter blue) band for each fit are 
included. Experimental data were taken from CAPRICE 31j, HEAT [30], 
AMS [33, MASS [3^, PAMELA [H] and ATIC-2 [35] M 

A.l Schematic graphics E' vs E. It shows how to determine integration bounds 
Cmaxjinin for a fixcd energy e' (in green). Maximum and minimum energies 
for the prime frame in terms of original frame energy (in red) are also shown. IllC 

C.l Schematic plot to illustrate how the charged pions decay produces polarized 
muons. Blue arrow represent left-handed polarization and red arrows right- 
handed polarized particles 111? 

C.2 Schematic plot x' versus 7. Integral limits on muon's 7 and frontiers between 

zones (red line), where spectrum changes behavior, are shown ll2C 

E.l Schematic view of creation of MoA~rule and how data splitting is performed. ll3C 
E.2 Example of rule created with MoA and the composition of the final rule by 
using linear interpolation. In this figure: red dots are the averages and red 
crosses correspond to extreme values. Note that final rule presents a non 

homogeneous distribution Il31 

E.3 Histograms made with a sample of 10^ events and 60 bins. An analytical 
distribution, MoA and uniform log-bin rule are plotted. The analytical dis- 
tribution is composed by three known distributions, a power-law (oc x"^) 
with two gaussian located in 10 and 100, all distribution have same weight. 
Binning rule obtained from MoA produce a more accurate solution than the 
uniform log rule, especially in the regions close to gaussians Il32 
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2.1 Parameters for charged Pion and Kaons inclusive CS according to Badhwar- 
Stephens-Golden parameterization [48|, I49j . The parameter units are such 
that the inclusive CS are in units of mb GeV~^ c'^ (Equations 12.41 and 12 . 10|) . 

2.2 Exclusive reactions in p-p collisions and the minimum configuration mass 
mx fEquation l2.9|l for charged pions and kaons [50 

2.3 Parameter for charged Pion and Kaons inclusive CS. These sets are for Tan- 
Ng parameterization [511I50| . The parameter units are such that the inclusive 
CS are in units of mb GeV^^ c'^ ( Equation 12. 11 1) 

2.4 Kinematics limits for Kamae et al. parameterization for electron and positron 
production 

2.5 Kamae et al. parameters describing electron spectra for arbitrary proton energy. 
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4.2 Typical combinations of cosmic ray propagation parameters that are com- 
patible with the B/C analysis [M]- The model MED has been borrowed 
from Table [37l1 Models Ml and M2 respectively maximize and minimize the 
positron flux for energies above 10 GeV - the precise extent of which de- 
pends on the mass of the DM particle, on the annihilation channel and also 
on the DM profile. Note that Ml is the same as MAX in Table \3~T\ but this 
is coincidental |65 

5.1 Values of parameters for Donato et al. and Shikaze et al. parameterizations 
of the proton and alpha particles fluxes. The units of fluxes for each param- 
eterization are m^^ sec^^ sr^^ (GeV/n)^^ 84 

5.2 Parameter values for the flt of the fc function [87 

5.3 Parameters values for best fits for electron flux of AMS fTT , HEAT [W and 
BETS ^33] experiments m 

D. l Most important charge kaon modes, probability and number of electron as- 

sociated to each of these [2] Il24 

E. l Number of bins estimation formulas for Sturgcs, Scott and Freedman- 

Diaconis models [1131 11141 1115j . For all models, n is the number of data, 
s^max / mill ^rc extreme values present, a is the standard deviation and IQR(a;) 
is the interquartile range function |12S 
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